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ABSTRACT 

SN 1987 A provides a unique opportunity to study the evolution of a supernova from explosion into very late phases. Owing to the 
rich chemical structure, the multitude of physical processes involved and extensive radiative transfer effects, detailed modeling is 
needed to interpret the emission from this and other supernovae. In this paper, we analyze the late-time (about eight years) Hubble 
Space Telescope spectrum of the SN 1987A ejecta, where 44 Ti is the dominant power source. Based on an explosion model for a 
19 M Q progenitor, we compute a model spectrum by calculating the degradation of positrons and gamma-rays from the radioactive 
decays, solving the equations governing temperature, ionization balance and NLTE level populations, and treating the radiative transfer 
with a Monte Carlo technique. We obtain a UV/optical/NIR model spectrum that reproduces most of the lines in the observed spectrum 
with good accuracy. We find non-local radiative transfer in atomic lines to be an important process also at this late stage of the 
supernova, with ~30% of the emerging flux in the optical and NIR coming from scattering/fluorescence. We investigate the question 
of where the positrons deposit their energy, and favor the scenario where they are locally trapped in the Fe/He clumps by a magnetic 
field. Energy deposition into these largely neutral Fe/He clumps makes Fe I lines prominent in the emerging spectrum. With the best 
available estimates for the dust extinction, we determine the amount of 44 Ti produced in the explosion to be 1.5^ j x 1CT 4 M Q . 

Key words. Line:formation - Line:identification - Radiative transfer - Supernovae:individual:1987A 



1. Introduction 

The distance of only ~50 kpc to SN 1987 A makes a number of 
unique observations possible. In particular, the supernova can 
after more than 20 years still be observed as resolved ejecta, 
glowing from radioactive input. About four years after explod- 
ing, the supernova entered the 44 Ti-dominated phase. Since then, 
its emission has been that of a cool (~10 2 K) gas powered by 
radioactive decays, complemented by freeze-out emission from 
the hydrogen envelope. Owing to the long life-time of 44 Ti (~85 
years), the character of the spectrum has changed little since this 
transition. 

To understand and interpret the late-time emission from the 
ejecta, detailed spectral synthesis modeling is needed. The spec- 
tral formation process in supernovae is complex, especially at 
late times when all important processes are non-thermal, but 
radiative transfer effects may still be important. High-energy 
gamma-rays and positrons from the radioactive decays deposit 
their energy into free and bound electrons, producing a popula- 
tion of fast electrons which heat, excite and ionize the gas. The 
ionizations produce secondary electrons that additionally con- 
tribute to these processes. The degradation process is quantified 
by solving the Boltzm ann equation for each of the zones present 
in the supernova (e.g. lXu & McCravlll991l; iKozma & Franssonl 
1 19921 KF92 hereafter). Once the deposition in the various chan- 
nels have been determined, the temperature, ionization structure 
and NLTE level populations can be computed by solving the 
equations for thermal a nd statistical equilibrium or their time- 
depend ent variants (e.g. lde Kool et alJl998tlKozma & Franssonl 
Il998allbl KF98 a, b hereafter). Because the ionization balance 
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affects the solution for the radioactive deposition, the equation 
systems are coupled and a solution is found by iteration. 

The equilibrium also depends on the radiation field, which 
should simultaneously be solved for. The strong velocity gradi- 
ents in supernovae imply unique radiative transfer effects. As 
they travel through the ejecta, photons are continuously red- 
shifting with respect to the comoving frame, and are there- 
fore exposed to absorption in lines over a broad wavelength 
range. Because especially iron-group elements have a large 
number of lines at UV/blue wavelengths, radiation here ex- 
periences a complex transfer and a significant fraction of it 
emerges by fluorescence as a quasi-continuum at longer wave- 
lengths dLi & McCravlll99"rj|) . The efficient blocking by atomic 
lines is equivalent to an effective continuum opacity known 
as line blocking. The emerging UV emission may often be 
due to 'holes' in this li ne blocking where radiation can escape 
dMazzali & LucvHl993l) . Although the optical depths decrease 
over time, the multi-line transfer is important for modeling su- 
pernova spectra for many years or even decades after explosion. 

The observational determination of the masses of the three 
main radioactive isotopes 56 Ni, 57 Ni, and 44 Ti constitutes one of 
the main constraints on explosion models of core collapse super- 
novae. The production of these isotopes is sensitive to both den- 
sity and temperature, and thus to the explosion dynamics (e.g. 
IWooslev & Hoffman! 1 99 ll) . During the first ~500 days the decay 
of 56 Co (which is the product of rapid 56 Ni decay) dominated the 
bolometric light curve of SN 1987 A. From the li ght curves, the 
56 Ni mass could be determined to be 0.069 M Q dBouchet et al.l 
119911) . The determination of the 57 Ni mass was complicated 
by the fact that when this isotope became important, the bolo- 
metric light curve was already affected by the delayed release 
of ionization energy (freeze-out). A time-dependent modeling 
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that took this into account resulted in an esti mated 57 Ni mass 
of -3.3 x 1(T 3 M Q dFransson & Kozmalll993l FK93 hereafter), 
which agrees well with the mass derived from observations of 
the inf rared Co H/Fe II lines and the gamma- rays emitted in the 
decay dVarani et al.lll990HKurfess et al.lll992l) . 

The determination of the 44 Ti-mass is complicated by a 
number of factors, the first being the freeze-out contribution 
just as for 57 Ni. Second, at late times the low temperature 
of the gas and the presence of dust lead to most of the de- 
posited energy being re -radiated in the largely unobservable 
far-infrared, which prohibits a determination of the bolometric 
luminosity. The 44 Ti-mass must therefore be determined from 
detailed modeling of the fraction of the energy that emerges 
in the (observable) UV/optical/NIR bands, and must consider 
the time-dependent freez e-out. From this type of modeling, 
iFransson & Kozmal d2002l FK02 hereafter) determined the 44 Ti- 
mass to (0.5 - 2) x 10" 4 M , which agrees with the ( 1 - 2) x 10 
M q estimated from a ne bular analysis of the eight-year spectrum 
by IChugai et al.l (1 19971 C97 hereafter). ILundqvist et all (1200 ll) 
derived an upper limit of 1.1 x 10 -4 M Q , based on the non- 
detection of the [Fe II] 26 //m line by ISO, which rested on the 
assumption of local positron deposition and the absence of any 
dust cooling. 

Another possibility for determining the 44 Ti-mass is by direct 
detection of the gamma-rays produced in the decay. Only in one 
case, Cas A, has such a detection been made dlvudin et al.ll 19941 
iVink et al.ll200lh. w ith the mass estimated to x 10" 4 M Q 

( Renaud et alJl2006h . For Cas A we do not know the masses of 
the 56 Ni and 57 Ni isotopes though. For SN 1987 A, the detection 
limits by INTEGRAL result ed in only a weak upper limit of 1 X 
10 3 M Q for the 44 Ti-mass dShtvkovskiv et al.ll2005h . 

In this paper we present a detailed model of the 
UV/optical/NIR spectrum of SN 1987 A at an age of eight years, 
using a 19 M explosion model as input. Using a new code 
with a detailed radiative transfer treatment, our objectives are 
to see if the spectral model can reproduce the main features in 
the observed spectrum, to understand the contributions by var- 
ious zones and elements to the spectral formation process, to 
quantify the importance of the (non-local) line transfer, and to 
determine the 44 Ti-mass as accurately as possible. Previous dis- 
cussions of the la te-time ejecta spectra can be found, e.g., in 
IWang et all d 19961) and C97. Our paper also serves to describe 
the code we developed for reference in future modeling. 



2. Modeling 

In KF98 a, b, a self-consistent model for the nebular phase spec- 
trum of SN 1987 A was presented. This included a detailed cal- 
culation of the gamma-ray/positron thermalization and a deter- 
mination of the time-dependent ionization, excitation, and tem- 
perature structure of the ejecta. Comparing the resulting spectra 
from two different explosion models, the evolution of all major 
emission lines were analyzed. IFransson & Kozmal d2002l) con- 
tains an update of this work, where the modeling of some of the 
most important lines and the broad band photometry were dis- 
cussed. 

A shortcoming in these models was that (non-local) radia- 
tive transfer in atomic lines was not included, which introduced 
an uncertainty for the internal radiation field and thereby for the 
ionization balance, and also for the emerging spectrum. Here, we 
calculate in detail the line scattering and fluorescence throughout 
the spectrum (Sect. 12.5b . which significantly improves the accu- 
racy of the model. In addition, we updated the atomic physics 



by adding new atomic data and including more and larger model 
atoms. 

While we do a more accurate radiative transfer calculation, 
this comes at the expense of time-dependent effects as our model 
assumes steady-state. In FK93 it was shown that in the envelope 
the recombination time scale became comparable to the radioac- 
tive decay time scale at ~800 days. From that time on, much of 
the ionization energy was not radiated instantaneously but on a 
longer time scale, and a time-dependent calculation was there- 
fore necessary. In the 44 Ti-dominated p hase the radioactive time 
scale becomes much longer, ~85 years dAhmad et al.l l2006T). and 
is less relevant. Instead, it is the adiabatic expansion time scale, 
i.e. the age of the supernova, that is the most relevant, and should 
be compared to the recombination and cooling time scales. In 
zones where a time-dependent treatment is necessary, we use the 
solutions obtained by an updated version of the time-dependent 
code used in KF98 a, b for the ionization balance. As we show 
below (Sect. 14.2b . this is only necessary in the hydrogen enve- 
lope. 

2.1. Ejecta model 

We base our ejecta model on the 1 9 M explosion model com- 
puted bv IWooslev & Hegerl (120071 WH07 hereafter). From the 
distribution of elements in this model, we construct seven types 
of zones; Fe/He, Si/S, O/Si/S, O/Ne/Mg, O/C, He, and H, where 
the zones are named after their dominant elements. The mass 
cuts for the zones are taken to be where the most common or 
second most common element changes. The mass of the Fe/He 
zone is adjusted to give a total 56 Ni mass of 0.07 2 M i n the 
ejecta (the value of 0.069 M from iBouchet etall d 19911) cor- 
rected for the slighly larger distance we use here), and the mass 
fractions of 57 Co and 44 Ti in this zone are then adjusted to give 
a total 57 Co mass of of 3.3 x 10" 3 (FK93) and a total 44 Ti mass 
of 1.5 x 10~ 4 M , which we later show to be our favored value 
for this isotope. The H zone mass is adjusted to match the ob- 
served estimate of ~8 M (KF98 b). Finally, we replace the He/C 
zone in the explosion model with a He/N zone, that has not ex- 
perienced any outward mixing of carbon. This is also based on 
modeling in KF98 b, where it was found that a significant He/C 
zone would result in [C I] lines much stronger than the observed 
ones. 

The masses and compositions of the zones in the model 
are given in Table Q] Abundances are taken as the value in 
the middle of each zone. While this does not exactly conserve 
the total element masses, the differences are small compared 
to using zone-averaged values, and the middle values should 
give a more consistent description for the typical composition. 
Because the WH07 models are based on solar metallicity pro- 
genitors, we replace the abundances of all elements with mass 
fractions differing by less than a factor 1.2 5 from the solar 
values with the LMC abundances taken from lRussell & Dopital 
(11992). A comparison of m odels with different metallicities in 
IWooslev & Weaver! dl995l WW95 hereafter) shows that most 
nucleosynthetic yields only weakly depend on metallicity, so the 
abundances for the processed elements should be very similar 
for solar and LMC metallicity progenitors. Observations of the 
circumstellar material indicate that CNO-burning products have 
been mixed out from the helium core into the hydrogen enve- 
lope dFransson et alJ 1989b . W e use the abundances of He, N and 
O for the envelope derived by IMattila et all (1201 Oh and of C by 
ILundqvist & Franssonl d 19961) (see TableQ}. 

We model the supernova as a spherically symmetric, homol- 
ogously expanding nebula, divided into a core and an envelope. 
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The core contains the heavy element zones (Fe/He, Si/S, O/Si/S, 
O/Ne/Mg, O/C), and fractions of the He and H zones, which 
are mixed inwards during the explosion. These fractions were 
chosen to be 25% (see Sect. 12.1. 11 1. The envelope contains the 
remaining 75% of the He and H zones. 

The extent of the core can be estimated from the widths of 
the lines that are emitted mainly from the heavy element zones, 
for example the iron lines and [O I] AA63Q0, 6364. The oxy- 
gen lines have maximum expansion velocities of ~ 1700 km s -1 , 
whereas the iron, cobalt, and nickel lines extend to ~2500 
kms" 1 , with s ome contribution even out to ~ 3500 km s -1 (see 
lMcCravlll993l for a review). Our single-core model forces a 
compromise value, which we choose as V core = 2000 kms -1 , 
the same as in KF98 a, b. A higher core velocity means less 
gamma-ray deposition, but because the positrons dominate the 
energy input at late times, this choice is not critical for the en- 
ergy budget here. 

2.1 .1 . The core 

We assume that the zones in the core experience complete 
macroscopic mixing, where they fragment into N c \ clumps that 
become completely mixed in velocity space (see Fig.Q]). Support 
for this m acroscopic mixing co mes from theoretical considera- 
tions (ejjjMu^lterjetj^ 

(e.g. [Herant & Benzl|l99l[ iMueUer et al.lll99lt iKifonidis et alJ 
I2006 ; Hamme r et al.l 2010 ) , the early emergence of X-ray s and 
gamma-rays dDotani etafl 119871; iPinto & Wooslevl 119881) . and 
from the similar line profiles of many of the elements (e.g. 
ISpvromilio et al. I fl990h . We use a va lue of jV c i — 1 00 based 
on the analysis of the iron clumps in iLi et al.l d!993l) . In Sect. 
14.61 we find that the model spectrum is insensitive to this 
choi ce. We further assume that no microscopic mi xing occurs 
(e.g. lFransson & Chevalierll989UFrvxell et al.lll991l) . so that the 
zones retain their compositional integrity. 

The clumpy structure of the core calls for a special handling 
of the computational grid. We use a gridless domain for the core 
region; when a clump emits photons, the position of the clump 
within the core is selected by a random draw. When we follow 
photons through the core (Sect. l2~5T l, exiting one clump leads to a 
random draw that determines which type of clump to encounter 
next. We set the probability of encountering any given type of 
zone to be proportional to the total surface area of that zone type. 
A random number also selects the impact parameter on the new 
clump. The path traveled through zone ; is then proportional to 
A, • Ri oc Vj, which recovers the desired property that a photon 
should spend its time in a given zone in proportion to the filling 
factor of the zone, e, . 

The core structure is a 'virtual' one in the sense that there is 
no fixed density grid; we have a set of spherical clumps which to- 
gether make up the total volume of the core, but which obviously 
cannot be exactly fit together since they are all spheres. There 
are several advantages to this virtual grid. The extension of the 
Monte Carlo philosophy to encompass also a randomization of 
the grid is a powerful and appealing way to make computations 
in clumpy media such as this. While the radiative transfer can 
still be done realistically by the use of surface and impact proba- 
bilities as described above, a virtual grid naturally offers a way to 
capture the situation of macroscopic mixing and to parametrize 
the degree of fragmentation of the core. Finally, the line profiles 
automatically become smooth rather than jagged, as when using 
shells. 

The filling factors for the oxygen zones are set from the 
oxygen number density derived from the evolution of the [O I] 



6300 A/6364 A line ratio n Q = 6.2 x 10 10 cm ' 3 at 100 days 
(ISpvromilio & Pintolll99U ILi & McCravlli992h . This allocates 
~10% of the core volume to oxygen clumps. 

The iron clumps are believ ed to expand because of their 
content of radioactive materials (IWooslevlll988HHerant & Benzl 
[T^lBaskol[T99l . which was verified for SN 1987 A by mod- 
eling of the cooling lines from the iron clumps dLi et al.lll993l) . 
The iron zone was found to have a filling factor of fFe/He <; 0.3, 
assuming Vc 0re = 2500 kms" 1 . The same density limit using 
Vcore = 2000 km s -1 corresponds to ep e /He <: 0.6. We use a value 
of 0.6 in our modeling here. 

In the explosion model, the Si/S zone contains ~5% of the 
56 Ni, and can therefore be expected to expand as well. For this 
zone, we therefore use a density in between the oxygen zones 
and the iron zone, chosen to be ten times lower than the oxygen 
zone density. This gives a filling factor of ~3%. In Sect. 14.3.21 
we find support from our modeling for a low density i n this zone. 

Hydrodynamical simulations by iFrvxell et al.l (Il99ll) find 
~1 Mo of hydroge n mixed inside 2000 kms 1 , while 
iHerant & Beiizl d 19921) find ~2 M . From modeling of the Ha 
line in SN 1987 A, KF98 b estimate the presence of 2.2 M of H- 
zone gas inside 2000 km s -1 . Here, we take 2 M Q of H-zone ma- 
terial to be mixed inside the core, which corresponds to ~25% of 
the total H zone mass. We assume the same fraction (25%) of the 
He zone to be mixed into the core as well, giving a He-zone core 
component of 0.3 M and a He-envelope component of 1.0 M Q . 
The total mass of the core, including the mixed-in He and H, is 
then 5.2 M Q . 

The filling factors for the H and He zones in the core are set 
by allocating the remaining volume so that these zones obtain 
equal number densities. This gives the H component a filling 
factor of ~26% and the He component a filling factor of ~2%. 
See also KF98 a for further discussions regarding filling factors. 



2.1 .2. The envelope 

Outside the core we attach an envelope consisting of an inner 
He zone followed by logarithmically spaced (constant dR/R) H- 
zone shells. As mentioned above, the total He-zone mass is 1.3 
M Q , of which we put 1.0 M in the envelope and 0.3 M Q in the 
core. The total H-zone mass is taken as 8 M from Ha modeling 
(KF98 b), of which we put 6 M in the envelope and 2 M in the 
core (see Sect. 12.1. Ti l. 

We use a density distributi on of the H-envelope based 
on hydrodynamical modeling by Shigevama & Nomotol d 19901) 
(model 14E1). The validity of this model is confirmed by mod- 
eling of the Mg II ,12800 and Mg 1 ,12852 features in C97. The 
density profile can be described as gradually steepening from an 
initial -2 power law to an asymptotic ~-8 for V > 5000 kms -1 . 
We terminate the envelope at a velocity of 10,000 km s -1 . 

2.2. Energy deposition and degradation 

We calculate the deposition of gamma-rays emitted by the 
clumps cont aining 56 Co, 57 Co and 44 Ti. Based on the sim- 



ulations by IColgate et all dl980l) : IWooslev et all dl989l) and 
iFransson & Chevalier! ( 19891) we assume absorption with ef- 
fective opacities of k( 56 Co) = 0.030 (f/0.5), k( 57 Co) = 
0.072 (I /0.5) and *( 44 Ti) = 0.040 (f /0.5) cm 2 g -1 , where Z and 

A are the average nuclear charges and atomic weights in each 
zone. We do not include other radioactive isotopes such as Co 
and 22 Na, whose abundances are highly uncertain. We discuss 
the possible effects of these in Sect. [5] 
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Table 1. Chemical compositions (mass fractions) of the zones used in the model. 



Zone 


Fe/He 




Si/S 




O/Si/S 




O/Ne/Mg 


O/C 




He 






H 






Mass (M Q ): 


0.092 




0.14 




0.16 




1.9 




0.58 




1.3 






8.0 






Mass fractions: 


56 Ni 


0.75 




0.022 

































57 Co 


0.034 




9.0 x 10- 


4 


1.0 x 10 


-5 


























44 Ti 


1.4 x 10 


-3 


5.2 x 10 


-5 


1.6 x 10 


-5 


























H 


2.3 x 10' 


-6 


4.2 x 10" 


8 


8.7 x 10 


■9 


4.6 x 10 


10 


2.1 x 10" 


-10 


2.5 x 


10 


-7 


0.60 






He 


0.14 




8.3 x 10~ 


-6 


4.3 x 10 


-(. 


2.0 x 10 


-6 


0.078 




0.99 






0.40 






C 


3.6 x 10 


-7 


1.4 x 10 


-6 


8.5 x 10 


-5 


0.031 




0.28 




3.0 x 


10 


-4 


2.4 x 


10" 




N 


1.4 x 10 


-6 


1.1 x 10 


-9 


1.8 x 10 


-5 


3.2 x 10 


-5 


4.7 x 10 


-6 


9.0 x 


10 


-3 


2.2 x 


10" 


-3 





9.7 x 10 


-(. 


1.6 x 10 


-5 


0.46 




0.70 




0.62 




1.6 x 


10 


-4 


1.8 x 


10" 


-3 


Nc 


1.0 x 10 


-5 


1.6 x 10 


-h 


1.4 x 10 


-4 


0.22 




0.018 




1.1 X 


10 


-3 


6.1 x 


io- 


-4* 


Na 


6.0 x 10 


-7 


1.0 x 10 


■6 


1.2 x 10 


-6 


7.0 x 10 


-3 


1.9 x 10 


-4 


1.8 x 


10 


-4 


1.6 x 


10" 


-5* 


Mg 


1.2 x 10 


-5 


1.6 x 10 


-4 


7.4 x 10 


-4 


0.042 




4.9 x 10 


-3 


5.4 x 


10 


-4* 


5.4 x 


10" 


-4* 


Al 


2.1 x 10 


-5 


2.6 x 10 


-4 


1.9 x 10 


-4 


4.0 x 10 


-3 


3.0 x 10 


-5* 


3.0 x 


10 


-5* 


3.0 x 


10" 


-5* 


Si 


1.6 x 10 


-4 


0.44 




0.27 




3.1 x 10 


-3 


2.8 x 10 


-4 * 


2.8 x 


10 


-4* 


2.8 x 


io- 


-4* 


S 


1.2 x 10 


-4 


0.39 




0.21 




2.6 x 10 


-4 


1.2 x 10 


-4 


1.2 x 


10 


-4* 


1.2 x 


10" 


-4* 


Ar 


1.2 x 10 


-4 


0.051 




0.043 




8.4 x 10 


-5 


5.4 x 10 


-5* 


5.4 x 


10 


-5* 


5.4 x 


10" 


-5* 


Ca 


1.3 x 10 


-3 


0.031 




0.012 




3.1 x 10 


-5 


2.3 x 10 


-5* 


2.3 x 


10 


-5* 


2.3 x 


10" 


-5* 


Sc 


2.3 x 10 


-7 


2.3 x 10" 


-7 


1.2 x 10 


-6 


1.6 x 10 


-(> 


5.6 x 10 


-7 


1.4 x 


10 


-8* 


1.4 x 


10" 


-8* 


Ti (stable) 


8.6 x 10 


-4 


2.9 x 10 


-4 


8.4 x 10 


-5 


6.2 x 10 


-6 


5.3 x 10 


-6 


2.3 x 


10 


-6 :|: 


2.3 x 


10" 


-6* 


V 


2.1 x 10 


-5 


1.3 x 10 


-4 


3.0 x 10 


-6 


6.4 x 10 


-7 


4.5 x 10 


-7* 


4.5 x 


10 


-7* 


4.5 x 


10" 


-7* 


Cr 


1.3 x 10 


-3 


3.2 x 10 


-3 


4.1 x 10 


-5 


1.3 x 10 


-5 


1.2 x 10 


-5* 


1.2 x 


10 


-5* 


1.2 x 


10" 


-5* 


Mn 


1.0 x 10 


-5 


6.6 x 10 


-4 


1.9 x 10 


-h 


3.8 x 10 


-6 


6.6 x 10 


-6 * 


6.6 x 


10 


-6 


6.6 x 


10" 


-6 * 


Fe (stable) 


2.3 x 10 


-3 


0.050 




4.3 x 10 


-4 


6.9 x 10 


-4 


7.2 x 10 


-4* 


7.2 x 


10 


-4* 


7.2 x 


10" 


-4* 


Co (stable) 


3.2 x 10 


-8 


5.5 x 10 


-9 


1.4 x 10 


-h 


1.9 x 10 


-4 


1.3 x 10 


-4 


2.4 x 


10 


-6 :\: 


2.4 x 


10" 


-6 :j: 


Ni (stable) 


0.037 




2.4 x 10 


-3 


1.2 x 10- 


-3 


5.6 x 10 


-4 


2.3 x 10 


-4 


4.8 x 


10 


-5* 


4.8 x 


io- 


-5* 


Mean atomic weight 


19.4 




31.6 




21.5 




16.9 




12.1 




4.0 






1.4 







Notes. Primordial LMC abundances are marked with an asterisk, and CNO abundances derived from observations of the inner circumstellar ring 
are marked with two asterisks. See also Sect. 12. II 



Table 2. Properties of the zones in the model. 



Zone 


Mass 


v m 


v out 


Filling factor (e) 


Density 


Rec. time 




(Mo) 


(kms" 1 ) 


(km s" 1 ) 




(cm" 3 ) 


(years) 


Fe/He 


0.092 





2000 


0.6 


1.8 x 10 4 


2.5 


Si/S 


0.14 





2000 


2.9 x 10~ 2 


3.5 x 10 5 


0.8 


O/Si/S 


0.16 





2000 


4.1 x 10~ 3 


4.2 x 10 6 


0.4 


O/Ne/Mg 


1.9 





2000 


7.3 x 10~ 2 


3.5 x 10 6 


0.9 


O/C 


0.58 





2000 


2.0 x 10~ 2 


5.6 x 10 6 


0.9 


He-core 


0.32 





2000 


1.5 x 10~ 2 


1.3 x 10 7 


0.5 


H-core 


2 





2000 


0.26 


1.3 x 10 7 


2.1 


He-env 


1.0 


2000 


2200 




1.7 x 10" 


1.3 


H-env- 1 


1.9 


2200 


2730 




2.5 x 10 6 


6.9 


H-env-2 


2.0 


2730 


3390 




1.4 x 10" 


10 


H-env-3 


1.2 


3390 


4210 




4.3 x 10 5 


17 


H-env-4 


0.61 


4210 


5230 




1.2 x 10 5 


29 


H-env-5 


0.24 


5230 


6490 




2.4 x 10 4 


49 


H-env-6 


0.083 


6490 


8050 




4.3 x 10 3 


105 


H-env-7 


0.028 


8050 


10000 




7.6 x 10 2 


332 



Notes. The last column shows the recombination time obtained from a steady state solution (see Sect. 14. 2b . 



From each zone containing radioactive material, we emit and 
track 56 Co, 57 Co, and 44 Ti gamma-ray packets as described for 
the UV/optical/NIR radiation in Sect. 12.51 with the difference 
that these packets travel in straight paths until they reach the 
edge of the nebula, depositing a fraction exp (—KpJ) of their en- 
ergy for each traveled distance /, where p\ is the density of zone i. 
We use radioactive luminosities taken from KF98, with the 44 Ti 
values updated for a life-time of 85 years instead of 78 years. 



After -2000 days, 44 Ti completely dominates the deposi- 
tion. The degradation of the deposited gamma-ray and positron 
energy into heating, ionizations and excitations is obtained by 
solving the Spencer-Fano equation using the routine developed 
by KF92. For this routine, we updated the collis ional cross 
sections for Fe II jRamsbottom et"ail I2005L 120071) and Ca I 
(ISamson & Berringtonl200i ). Ionizations are generally assumed 
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Fig. 1. Schematic structure of the used ejecta model. The core 
consists of several chemically distinct zones distributed as N c i 
clumps each. The envelope consists of an initial He zone (gray) 
followed by shells of H-rich gas (white). The densities of the H 
shells are decreasing outward. 

to leave the ions in their ground states, except for O II, for which 
we compute specific rates to the first few excited states. 

In our standard model, the positrons are assumed to be ab- 
sorbed on the spot. As noted before, e.g. in C97, this requires a 
magnetic field. The average positron energy in the 44 Sc decay is 
0.6 MeV and the ma ximum energy of the positrons is 1 .47 Me V 
dBrowne et al. I fl986h . The energy loss per unit length, dE/dx, 
is dominated by excitations and ionizations and is given by the 
Bethe formula 



dE 
dx 



Anrlm e c 2 p Z 
m u /3 2 A 



(I) 



where ro is the electron radius, f3 = V/c, I is the effective ion- 
ization potential of the atom, f z is a relativistic correction factor, 
and all other symbols have their usual meaning. In Fig. [2] we 
show the stopping range, p£/^f (g cm -2 ), for the most impor- 
tant elements in the core together with the total range, averaged 
over the different zones in the core. 

These ranges should be compared to the mass column den- 
sity D of the various zones in the core and envelope. For the total 
core, this is given by 



D = 0.93 



'M a 



i5M, 



,j\ 2000 km s- 1 / \8yrsj 



gem 



-2 



(2) 



The column density for a core with V core = 1800 kms~'at 7.87 
years is shown as the upper dashed line in Fig. [2] For comparison 
we also show the column density of the Fe/He zone as the lower 
horizontal line. 

That the column density of the Fe/He zone is nearly two 
orders of magnitude lower than the stopping range for Fe and 
He means that in the absence of a sufficiently strong magnetic 
field the positrons will not be trapped in this zone, even if it is 
not strongly fragmented into many small pieces. The figure also 
shows, however, that the total core has a sufficiently high col- 
umn density to trap the positrons even without a magnetic field. 
There should therefore be no positrons entering the envelope at 
this or earlier epochs. 

The presence of magnetic fie lds in sup e rnova ejecta has been 
sugges ted many times, see e.g. lAxelrodl (Il980l) . IColgate et~aTl 
(1 1 980h argue that any magnetic field present at the time of the ex- 
plosion becomes radially 'comb ed' and will therefore not be im - 
portant for the trapping (see also lRuiz-Lapuente & Spruitll998l) . 




E (MeV) 

Fig. 2. Positron stopping range for different elements as a func- 
tion of kinetic energy. The solid black curve denoted 'tot' shows 
the range weighted by the different zones in the core. The lower 
horizontal dashed line marks the column density of the Fe/He 
zone (assuming it to exist as a single clump), and the upper that 
of the full core for V C0le = 1800 kms . The long-dashed curve 
shows the cumulative positron distribution. 



While the positrons should be confined to the core at eight years, 
it is therefore not clear whether they are absorbed by the zones 
producing them (which is mainly the Fe/He zone, but also the 
Si/S zone), or if they can propagate into the other core zones be- 
fore they are absorbed. We will later investigate this question by 
comparing the spectra produced in the different scenarios (Sect. 

2.3. Statistical and thermal equilibrium 

For each zone we solve the equations for statistical equilibrium 
with respect to the ionization and excitation structure, as well as 
the equation of thermal equilibrium. These equations are essen- 
tially the ones described in KF98 a, with some modifications and 
improved rate calculations described below. For each iteration, 
the algorithm deployed is to first compute the ionization balance, 
then iterate solutions to excitation structure and thermal equilib- 
rium until the temperature stabilizes, and finally to calculate the 
radiative transfer (Sect. 12.5) . From the radiative transfer, new 
photoionization rates, photoexcitation rates (although these are 
set to zero in this work, see Sect. I2.3.2I I. and radiative heating 
rates are calculated, and used for the new solutions for ioniza- 
tion, excitation and temperature in the next iteration. Iteration 
continues until the maximum change in any zone temperature, 
electron fraction or ion abundance is less than 1%. We then also 
verify that the change in the emerging radiation field is neglegi- 
ble. 

2.3.1. Ionization balance 

Steady-state for the ionization balance is appropriate as long as 
the recombination time-scale (l/n e a) is short compared to the 
evolutionary and radioactive time scales. For the densities at the 
epoch studied here, this is satisfied (Sect. I4.21 i for all core zones 
as well as for the envelope helium zone. The hydrogen enve- 
lope, however, experiences freeze-out after ~800 days (FK93). 
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For the envelope, we therefore use the solutions for the hydro- 
gen ionization fraction from a time-dependent calculation, based 
on the model in KF98 a, given in Table [3] We do, however, use 
steady-state solutions for the other elements in the envelope, be- 
cause the time-dependent model lacks photoionizations by line 
radiation, and may therefore significantly underestimate the ion- 
ization degree for some of them. We also take the hydrogen en- 
velope temperatures from the time-dependent solution, because 
temperature and ionization balance are coupled. 



Table 3. The degree of hydrogen ionization, and the tempera- 
ture in the hydrogen envelope at day 2875, taken from a time- 
dependent calculation based on the model in KF98 a, b. 



Velocity 


Number density 


n(H II)/n(H I) 


Temperature 


(kms" 1 ) 


(cm" 3 ) 




(K) 


core 


1.7 x 10' 


6.5 x 10-' 


119 


2000 


3.2 x 10 6 


2.4 x 10- 4 


114 


2400 


2.5 x 10 6 


2.5 x 10- 4 


101 


2800 


1.5 x 10" 


3.3 x 10- 4 


93 


3300 


5.4 x 10 5 


7.3 x 10- 4 


89 


4100 


2.4 x 10 5 


1.3 x 10- 3 


85 


4440 


1.2 x 10 5 


2.2 x 10- 3 


84 


5000 


4.6 x 10 4 


4.4 x 10- 3 


84 


5500 


2.2 x 10 4 


7.5 x 10- 3 


83 


6000 


1.1 x 10 4 


1.2 x 10- 2 


81 


6500 


5.4 x 10 3 


1.6 x 10- 2 


76 


7000 


3.0 x 10 3 


2.0 x 10- 2 


73 


7500 


1.6 x 10 3 


2.4 x 10- 2 


70 


8000 


9.5 x 10 2 


2.8 x 10- 2 


67 


8500 


5.7 x 10 2 


3.2 x 10- 2 


64 


9000 


3.5 x 10 2 


3.5 x 10- 2 


64 


9500 


2.2 x 10 2 


3.8 x 10- 2 


62 


10,000 


1.4 x 10 2 


4.2 x 1Q- 2 


60 



Notes. Note that the zones here are somewhat differently structured 
than those we use for the calculations in this paper (Table |2}. 



For all other zones, the ionization balance is computed from 
the steady-state variant of Eq. (10) (and its generalized version) 
in KF98 a. We include the first three ionization stages of all el- 
ements listed in Sect. 12.3.21 The main improvement in our treat- 
ment here is that we now compute detailed photoionization rates 
from the Monte Carlo simulation (Sect. l2~5T l. For zone i and atom 
k, the photoionization rate is 



Z;M,k,j 



(3) 



k-J 



where j is the level, ATyy is the total number of photoionizations 
for zone i, atom k, level j, in the Monte Carlo simulation, V; is 
the zone volume and n* k • is the previous solution to the number 
density of level j. 

2.3.2. Excitation structure 

Given the ion abundances, steady-state is appropriate for the 
atomic level populations since the collisional and radiative time- 
scales are short. We compute NLTE solutions for the neutral and 
singly ionized elements of H, He, C, N, O, Ne, Na, Mg, Al, Si, S, 
Ar, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, and Ni, as well as for doubly 
ionized Fe, using model atoms that include fine-structure lev- 
els. With their large number of resonance lines in the UV and 



optical, iron-group elements are especially important for the ra- 
diative transfer, which motivated us to include also those with 
small abundances. Appendix lAl gives references for the atomic 
data that we use. 

For the line opt ical depths, we use the Sobolev approxima- 
tion (ISobolevll 19571) 



with the corresponding escape probabilities fif. 



1 



(4) 



(5) 



where i and j are the lower and upper levels, respectively. The 
Sobolev approximation is valid as long as the region of line 
interaction is much smaller than the region over which phys- 
ical conditions change, which is generally fulfilled in super- 
novae. The errors caused by line-ove rlap should be small in 
the nebular phase dLi & McCravlll99"6l) . For lines with small 
we also take into account continuum de struction probabilities 
(iHummer & Rvbicki|[T985l IChugailll987h . so that the effective 
escape probability is 



(6) 



where the probabilities are calculated as in KF98 a. The ef- 
fective radiative deexcitation rates used in the NLTE equations 
are then Aj^ ff . 

The statistical equilibrium is obtained by solving Eq. (15) in 
KF98 a. Deexcitation of the He I(2 3 S) state also occurs through 
Penning ionizations of H I. Our treatment here differs mainly in 
that we now compute detailed photoionization rates Pyy from 
the Monte Carlo calculation (Sect. I2.5l l 



. M.k.j 

1 i,k,J 



(7) 



As in KF98, we include non-thermal ionizations and pho- 
toionizations only to the first level of the ionized element. The 
only exception is for O II, where non-thermal ionizations to the 
excited states 2 D , 2 P a and 4 P are also included. Non-thermal 
ionizations are included only from the ground multiplets in the 
atoms, whereas photoionizations can occur from excited states 
where cross sections are available. 

In this work, we set all photoexcitation/deexcitation rates to 
zero. Although these can in principle be computed from the ra- 
diative transfer calculation and be included in the population 
equations, this would force a slow single-stepping Monte Carlo 
transfer. We also aim to study the scattering/fluorescence process 
in detail (Sect. 14.3b . which is easier if we can track the Monte 
Carlo packets through all scattering events until they escape. 
Finally, this approximation also allows us to compute faster (ap- 
proximate) NLTE solutions including a fewer number of states. 
At these late times, photoexcitations do not strongly influence 
the solutions beyond the direct scattering/fluorescence. 

For the treatment of charge transfer, see Sect. 12.41 

2.3.3. Thermal balance 

For the core zones, the temperature is computed by balancing 
heating with cooling. Steady-state is appropriate as long as adia- 
batic cooling is neglegible. This may not be fulfilled at late times 
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for the oxygen, helium, and hydrogen zones (KF98 a). However, 
we compared the steady-state temperatures in these zones with 
those obtained from the time-dependent model (KF98 a, b), and 
they are found to be very similar. We therefore use the steady - 
state approximation for all core zones. For the envelope zones, 
we use the temperatures calculated in the time-dependent model. 
These are given in Table [3] 

Heating has contributions from non-thermal deposition, pho- 
toionizations, collisional deexcitations (occuring in the Monte 
Carlo calculation), charge transfer reactions, free-free absorp- 
tions, and Penning ionizations 



H - H nt + H v \ + // c d + H ct + fig + H pe 



(8) 



The terms // p ;, H c a and Hg are computed from the Monte 
Carlo simulation (Sect. 12.5b by summing all contributions 
from Eqs. ( l23l l. (l24l and (1341 . The charge transfer heating is 
dKingdon & Ferland1l996h 



H ct = ^«k«k<^jA£ kik < 



(9) 



k.k' 



where is the reaction rate between atoms k and k! ', and AE^> 
is the energy defect of the reaction, which is negative for en- 
dothermic reactions, which act as coolers. 

T he hea ting from Penning ionizations is, based on the rate in 
iBeU (fl970h 



//pe = 7.5 X 10 21 n He I(2 3 S) "HI • 



(10) 



Cooling is computed from Eqs. (3)-(6) in KF98 a. It has con- 
tributions from line cooling, recombination cooling, and free- 
free emission. 



2.4. Charge transfer 

We u s e charge transfe r rates taken from Rutherford etahl 



d 19711). Rutherford et al. 


( 


1972 


). lArnaud & Rothenflud (11985). 


Pequianot & Aldrovandi 


( 


1986 


), iKimura et all dl993l). Swartzl 


(1994), 


Kinsdon & Ferlandl (1996), Stancil et al. (1998t) and 


Zhao et al.l (12004). For reactions where no measured or calcu- 



lated rates exist, we use the rec ipe for estimating rates given by 
iPequignot & Aldrovandi d!986l) . 

The many unknown or uncertain rates introduce an uncer- 
tainty for the level of ionization for some of the elements. For 
some zones, this prohibits a definite determination of which el- 
ements that re-emit the non-thermal and radiative ionization en- 
ergy. We discuss this and the potential effects on the spectrum in 
Sect.l431andl4771 

When the final states are unknown, we assume these to be 
the ground states of the elements. An exception to this is the 
important OII + CI— > O I + C II reaction, where we assume that 
the oxygen atom is left in the excited 2p('D) state, which is the 
one closest to resonance with an energy defect of 0.4 eV (Sect. 
14.2.4b . We also assume that the reaction Ca I + O II -> Ca II + 
O I o ccurs to the excited 5p( 2 P) state in Ca II dRutherford et al.l 
119721) . 

2.5. Radiative transfer 



Radiative transfer is treated with a Monte Carlo technique. The 
application of this method for dealing with differ entially expand- 
ing as trophysical envelopes was described by I Abbott & Lucv 
1985 ), whe re it was applied t o the case of stellar winds. iLucv 



computing photospheric-phase supernova spectra. The main dif- 
ficulty working with the equation of radiative transfer is the 
complexity introduced by the line overlaps that are caused by 
the differential velocity field. This gives a non-local coupling 
for which a Monte Carlo treatment is better suited. In addition, 
Monte Carlo radiative transfer is easily generalized to 2D and 
3D, and is well suited for parallelization. 

The early Monte Carlo codes mentioned above allowed for 
electron and resonance line scattering, where deexc i tation s oc- 
cured in the sa me transition as the absorption. iLucvl (I 1999b and 
iMazzalil (|2000) improved on this treatment by including fluores- 
cence in_all_ttamitions directly connected to the (initial) upper 
level . iKasen et"al ] d2006h used the same treatment for the fluores- 
cence in computing photospheric-phase, time-dependent spectra 
for Type la supernovae. To allow for a more exact treatment of 
the fluorescence, but still retain the properties of indivisibility 
and indestructibility for the photon packets, a Monte Carlo for- 
malism based on th e concept of macro-atoms was developed in 
a series of papers bv ILucvl (120021 120031 120051) . Applications and 
further developments of this method are described in e.g. ISiml 
(l2007h . [Sim"et al.ld2007l) and lKromer & Simld2009l) . 

Here, we implemented a Monte Carlo algorithm with a sim- 
ilarly high level of physical realism by computing the complete 
deexcitation cascade following each photon packet absorption. 
A new photon packet is created for each radiative deexcitation 
in the cascade, with collisional deexcitations being allowed for 
as well. The (potential) creation of several packets following the 
absorption of one packet is handled by a recursive technique. 
We find that this works well and have therefore not enforced the 
property of indivisibil ity otherwis e commonly used to simplify 
the bookkeeping (e.g. lLucvll2002l) . Also, the property of inde- 
structibility devised in those papers has its main advantage as 
a Lambda accelerator in the photospheric phase, and is not en- 
forced here. 

2.5.1. Packet emission 

Below, quantities are denoted by primed symbols in the comov- 
ing frame, and by unprimed symbols in the rest frame (observer 
frame). We ignore terms of order (V/c) or higher in the trans- 
fer, except for the critical Doppler shifts. The letter z denotes a 
random number between and 1, newly sampled for each com- 
putation. 

Total emissivities j' A ,(i,j) from lines and continua (two- 
photon, free-bound and free-free) are obtained for each zone i 
and wavelength bin j from solutions to the ionization balance, 
temperature and NLTE level populations (Sect. 12.3b . Photon 
packets are then created with energy 



(ID 



19871) and iMazzali & Lucvl d 19931) describe the method for 



where V is the volume of the zone, A/lj is the width of the wave- 
length bin, and iVy is the number of packets to split the emis- 
sion into. The wavelength A' is chosen as the maximum wave- 
length of bin j. This choice ensures that no line self-absorptions, 
which are already included in the Sobolev escape formalism, oc- 
cur in the transfer. We currently use a uniform weighting for 
the number of packets emitted per zone and wavelength bin, 
M,j = Nq. From inspection of the emerging spectrum (which we 
also smooth, see Sect. I2.5.4i i, we find iVo ~100 to be sufficient 
for obtaining neglegible Monte Carlo noise, for binning width 
AA'/A' = 1Q- 3 . 
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The starting radius of the packet is sampled from 



or ■ 



(12) 



which corresponds to uniform sampling in a shell bounded by 
fi,in = V'i inf and ri jOUt = Vj, ou tf. The shell velocities are listed in 
Table ' 

Isotropic emission corresponds to sampling a direction co- 
sine fi from 



fi = -l+2z. 



(13) 



If the zone is a core zone, the internal radius in the (spherical) 
clump is also determined from 



r - 7 1/3 ff- , 

' int — "l, clump 



where /?i, c iump is the radius of the clumps of zone type i. An in- 
ternal direction cosine yU; nt is also randomly drawn according to 
Eq. O- 



2.5.2. Packet propagation 

For each packet to be emitted, a random number is drawn to 
determine at which optical depth the packet should be absorbed 
(if ataU) 



Tabs = -ln(l -Z) . 



(15) 



Next, the distance to the next critical point is computed. The 
comoving wavelength of the photon packet is continuously red- 
shifting as the packet travels through the differentially expanding 
ejecta. From the NLTE solutions, we have a wavelength ordered 
list of all lines with optical depth Ty > 10~ 3 . The velocity dis- 
tance to the next line is computed as 



AV lme = c(AljA'-l) , 



(16) 



where is the (atom rest frame) wavelength of the the line 
closest to A' in that list. 

The distance to the next ionization edge is computed as 



ion. edge 



= ^(<,edgeM'- 1 )' 



(17) 



where A , is the (atom rest frame) ionization edge closest in 

ion. edge v / o 

wavelength to X. Also, the distance to the next shell AV sn eii is 
computed, as is the distance to the clump edge AV c i ump , if the 
zone is a core zone. The velocity distance to the next critical 
point is then chosen as the smallest of these four distances 

AV = min{AVW, AV ion . edge , AV she ii, AV clump } . (18) 

The continuum optical depth for this distance is computed as 



= AVt 



ffi,dust + 2 a^(X) + a m {X) 



(19) 



where a^dust is the absorption coefficient for the dust 



Ku, < v C0K 

otherwise, 



T d 

oruus, - — 

y core 1 

= otherwise, (20) 

a^'\X) is the photoionization absorption coefficient, and 
&i,s(X) is the free-free absorption coefficient. Electron scatter- 
ing is neglegible at the epoch studied here, but can easily be 
included as well. 



The accumulated optical depth r tot is increased by r cont . 
If this causes T tot to transcend T aDS , the packet is continuum- 
absorbed, and propagation terminates. The total number of pho- 
toionizations for zone ;, atom k, level j, is then increased by 



AM,k,j = 



^"cont 



N, 



packet 



where Af pac ket is the number of photons in the packet 
E' 



packet 



he I A' 



(14) The photoelectric heating rate contribution is 



(21) 



(22) 



(23) 



where is the ionization potential for atom k, level j. The free- 
free heating rate contribution is 



1 AVta iS (X) , 
AH' S = — ^-Le' 

y'i T cont 



(24) 



If continuum absorption does not happen, the position, flight 
angle, wavelength and energy of the packet are now updated. 
The new position is 



1/2 



f-new = [r 2 + (AVt) 2 + 2rAVtfi] 
The new flight angle is 



sin y new = sin V . 

'"new 

The new comoving wavelength is 

^new = A' |l + — j , 

and the new comoving energy is 
E nf >w — E ( 1 + 

c 



(25) 



(26) 



(27) 



(28) 



If the new point r new corresponds to the passage of an ioniza- 
tion edge or a shell edge, new velocity distances are computed, 
starting from Eq. ([ToT l. and the procedure is repeated. If r new cor- 
responds to entering a new core clump, its type is chosen from 
surface impact probabilities given by 



R 2 , 

i, clump 
"7 ^j, clump 



(29) 



which is valid as long as all clumps have the same number of 
fragments {N c \). An impact parameter is also drawn as 



/^imp — Z 



(30) 



Propagation then continues into the new clump, starting from 
Eq. ©. 
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2.5.3. Line interaction 



^obs 



If ?"new corresponds to reaching a line with ry > 10 , T tot is 
increased by rii ne . If T tot then transcends r a b s , the packet is ab- 
sorbed in the line. 

All downward radiative and collisional transition rates are 
then computed. If u is the upper level of the line, the radiative 
deexcitation rates are (i < u) 

Rv=Avj0£. (31) 

The collisional deexcitation rates are (lOsterbrock & Ferlandl 
120061) 



_ 8.629 x 10~ 6 T u ,j(r) 
T 1/2 Xn 



(32) 



where T u ;(r) is the effective collision strength. The rates are 
normalized to obtain the relative transition probabilities, and one 
of the deexcitation channels is selected by a random draw. 

If this is a radiative transition (with a wavelength shorter than 
a pre-selected cut-off value, chosen as 2. 1 fim here) a new photon 
packet is emitted. The energy of this packet is 



-^out — ^new (A-Eui/^u) > 



(33) 



where A£ u ; is the energy difference between levels u and i and 
E u is the energy of level u. A direction is randomly drawn from 
Eq. dT3l ). and the packet is followed, starting from Eq. (Tl5t . until 
it (and its potential offspring packets) are absorbed or escape the 
ejecta. 

If the deexcitation is instead a collisional one, a contribution 
to the volumetric heating rate 



(34) 



is recorded. We ignore the possibility of collisional excitations 
in the cascade. This is a good approximation at late times when 
the low temperatures make upward collisions much rarer than 
downward ones. 

New deexcitation probabilities from level / are now com- 
puted from Eqs. (fJTJ and d32l i. and a channel is again randomly 
selected. The process repeats until the atom has deexcited to the 
ground state. The energy in each deexcitation from level i to level 
j is given by Eq. ( T33T > with AE U [ replaced by AE-,j. This treatment 
means that 



2> 



F' 



(35) 



While this becomes only an approximate treatment for the en- 
ergy distribution of the fluorescence when the absorption occurs 
from an excited state, this is rare at late times when most atoms 
are in their ground states. This treatment will likely be modified 
for calculations at earlier epochs when this may not be the case. 
Note that the outgoing energy (assuming only radiative deexci- 
tations) is generally not conserved in the rest frame, although it 
is in the comoving frame. This corresponds to adiabatic losses 
of the radiation field. 



2.5.4. Emerging spectrum 

When it reaches the edge of the nebula (with expansion velocity 
V e ), the observed wavelength and energy of the escaping packet 
are obtained as 



lobs 



r (.-**) 



(37) 



(36) 



where fi e is the direction cosine of the packet upon crossing the 
outer edge. The escaping packets are binned according to A b S , 
using a wavelength bin size of AA/A = 0.001. The final spectrum 
is obtained by smoothing this distribution with a box-car filter of 
width AA/A = 0.002. 

An alternative method for computing the emerging spectrum 
would be to solve for full-size atoms, including photoexcita- 
tion rates in the NLTE level population solutions, and compute 
the formal integral. However, computing large NLTE atoms is 
time-consuming and also numerically challenging in terms of 
convergence and accuracy. The line transfer mainly depends on 
the populations of the ground states and other low-lying states, 
which are more reliably computed than high-lying states dLucvl 
1 1 999b . We therefore choose to compute smaller NLTE solutions 
without photoexcitation/deexcitation rates and treat the scatter- 
ing/fluorescence process explicitly in the Monte Carlo transfer, 
which also allows us to investigate its properties from the emerg- 
ing packet histories. 

Because our model assumes steady state, we checked that the 
typical photon flight times f p are short compared to the dynamic 
time scale. Crossing the ejecta takes a time of about ~ V e t/c for 
the photons, which means that f p /f ~ V e /c <k 1. However, UV 
photons may experience multiple line scatterings before emerg- 
ing and their flight times could therefore become long. Here, the 
presence of dust in the ejecta limits the total flight times to a 
maximum of a few crossing times, and steady state is therefore 
justifiable from the radiative transfer perspective as well. 

2.6. Dust 

There is strong evidence from blue-shifting line peaks, a devel- 
oping IR excess, and a simultaneous decrease in the optical flux 
that dust formed in the inner parts of the ejecta, starting around 
day 530 dLucv et al.l 1 19891: IWooden et al.lfl993h . or possibly as 
early as day 350 dMeikle et al .1119931). By day 700 most of the 
dust formation was probably over dLucv et al.1119891 The appar- 
ently weak wavelength dependence of the dust extinction and 
the absence of any dust features in the IR emission suggest that 
the dust was at l east partly formed in optically thick clumps 
dLucv etal.ll 1991) . 

The last far-IR observations are from day 1731 
dBouchet et alj 11996 ), when the dust temperature was esti- 
mated to be ~ 155 K. Fassia et a fl (120021) find that essentially al I 
the NIR ejecta lines have blue-shifted peaks even at day 2112, 
which indicates that the dust was still optically thick at that time. 
We determined the median value o f their measured lin e-shifts to 
500 kms -1 . This agrees well with I Wang et al.l d!996l) . who find 
blue-shifts of ~400 kms -1 for several optical lines at 1862 and 
2210 days. 

Between 2000 and 6000 days there were no i nstruments ca- 
pable of detecting the FIR emission by the ejecta. iBouchet et al.l 
(120041) identify a Gemini N filter detection at 6094 days with 
ejecta dust. Observations at 6526 days were dominated by the 
heated dust emission fro m the ejecta - ring co llision, and no 
ejecta dust was detected dBouchet et alj 120061) . But the HST 
imaging shows a clear 'hole' in the centra l parts of the ejecta 
even at these late times dLarsson et al 1120111) . suggesting that the 
dust is still optically thick. This may, however, a lso be a result 
of ext ernal X-ray heating from the ring collisions dLarsson et al.l 
120111) . 
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We examined the observed line profiles of Mg I] /14571 and 
Ha in the eight-year spectrum (see Sect. |3]for a description of the 
observed spectrum). Both show distinct blue-shifts of their peaks 
of ~700 and ~600 km s _1 , respectively, indicating that the dust is 
still providing extinction at this epoch. This was also noted from 
the Ha profile by C97. The constancy of the line shifts compared 
to earlier epochs suggests that the dust resides in co-expanding 
clumps with optical depths much higher than unity also at this 
late epoch. 

Assuming emission and absorption in a homogeneous sphere 
of velocity Vd, the dust radial optical depth Td c an be found from 
the blue-shift of the peak A V dLucv et al.ll 1 989h 



AV ln(l+T d ) 

V d + T d 



(38) 



Using this equation for the Mg I] /14571 and Ha line shifts with 
Vd = 2000 km s 1 gives Td = 1.0 — 1.2, similar to it s value 
at earlier epochs dLucv et al.ll 19891 IWooden et"ai1ll993l) . Using 
AV = 400 km s" 1 , as f ound from the day 1862 and day 2210 
spectra by IWang et al.l ([l996), gives a val ue of Td _ 0.6, an d 
using AV = 500 kms -1 , as determined in [Fassia et al.l (|2002|) . 
gives Td = 0.8. 

Because it is unknown exactly in which zone(s) the dust re- 
sides, we model the dust as a constant, gray opacity within the 
core. From the estimates above, Td is likely to be in the 0.6 — 1.2 
range. For our standard model we choose Td = 1. 

We also note that dust in the ejecta can have strong effects 
on the cooling (KF98 a, b), with consequences for the spectrum 
depending on where the dust formed. 



3. The observed spectrum 

For the UV/optical range, we use the post-COSTAR HST Faint 
Object Spectrograph observations obtained on January 7 1995 
(2875 days or 7.87 years after explosion), which covers the 
wavelength range from 1650 to 9200 A. Details on the obser- 
vations and the data reduction can be found in C97. The flux 
calibration is stated to be good to ~5%. The aperture excludes 
the inner circumstellar ring, but the northern outer loop passes 
through the same line of sight as the central debris, and therefore 
contaminates the spectrum with some narrow lines. This is not a 
major problem because these lines can be easily identified. The 
strongest ones are C III] A 1909, [O II] A3 727, fO III] ,114959 , 
5007, [N II] ,U6548, 6584, and Ha ,16563 dPanagia et al.ll 19961) . 
We did not remove these narrow lines from the observed spectra 
displayed in the paper. 

For the near-i nfrared, we use the AAT spectrum from 
iFassia et"ai] d2002l) from day 2952. The IR spectrum contains 
contributions from both the ejecta and the inner equatorial ring. 
Because of the limited spectral resolution of A/AA = 120, nar- 
row lines from the ring can appear almost as broad as the ejecta 
lines. In addition, the removal of both continuum and line emis- 
sion from star 3 and possibly star 2 adds a significant uncertainty. 
The accuracy of the absolute flux levels is stated as +40%. The 
NIR spectrum is therefore useful for line identifications, but only 
for approximate quantitative comparisons. 

To compare th e model to the observations, we adopt a dis- 
tance of 51.4 kpc dPanag ia 2005), the extinction curve (for A < 
8000 A ) derived from the neighboring star 2 in IScuderi et al.l 
(119961) . which has E H v = 0.19. F o r A > 8000 A we use 
the extinction curve in Cardelli et al.l (fl989l A different anal- 
ysis bv lFitzpatrick & Walbornl (1 19901) find £ B -v = 0.16, and we 



take +0.03 to be a reasonable estimate for the extinction uncer- 
tainty. This dereddening gives a total luminosity in the 1600- 
8500 A range of 3.4 x 10 35 ergs" 1 . The value found by C97 is 
3.5 x 10 35 ergs -1 , using a distance of 50 kpc and £b-v = 0.2. 
We also adjust the model spectrum for interstellar UV absorption 
lines based on lBlades etall dl988l) . 

The dis t ance u ncertainty to SN 1987A is stated as 1.2 kpc 
in iPanagial d2005l) . which corresponds to a flux uncertainty of 
+4.8%. Dereddening the spectrum with £b-v = 0.16 or 0.22 
instead of 0.19 implies as change of ~ ± 13% in the total lumi- 
nosity in the 2000-8000 A range (which is the range we later use 
for estimating the 44 Ti mass). Combined, the observed luminos- 
ity then has an estimated uncertainty of ~ (5 2 + 4.8 2 + 13 2 )'^ 2 = 
+ 15%. 



4. Results 

4. 1 . Energy deposition and degradation 

At eight years, 44 Ti completely dominates the energy input over 
56 Co and 57 Co. Using a 44 Ti mass of 1.5 x 10~ 4 M , which we 
later show to be the best value for reproducing the observed spec- 
trum, the total energy deposition in the ejecta is 1 .9x 10 36 erg s" 1 . 
Of this, ~85% is positron energy; the other ~15% is gamma-ray 
energy, which represents ~5% of the 5.2 x 10 36 ergs^'gamma- 
rays emitted. The other ~95% of the gamma-rays escape the 
ejecta. 

The dominance of positrons, combined with the assumption 
that these are locally absorbed, means that most of the energy 
(~80%) is deposited in the Fe/He clumps. Most of the gamma- 
rays are absorbed in the H zones, giving them ~8% of the total 
energy deposition. The remaining ~12% are split between the 
silicon, oxygen, and helium zones. All energy depositions are 
given in Table 



4.2. Temperature, ionization, and emissivities 

The solutions for the temperature and ionization in the various 
zones are given in Table|5] From these steady-state solutions, we 
need to check the assumption that the recombination time-scale 
is shorter than the time-scale for density change njh — t/3 ~ 2.6 
years. The last column in Table [2] lists the recombination times 
( 1 / [n e a]) in years for the dominant element in the various zones. 
We see that a steady-state treatment of all core zones as well as 
of the helium envelope zone is justified, whereas the hydrogen 
envelope is experiencing freeze-out as expected. 

The temperature is 70-170 K in all core zones (Table|5]l. The 
temperatures could be even lower if molecular cooling or dust 
cooling are important. This would not have any strong influence 
on the UV/optical/NIR spectra, however, because non-thermal 
processes dominate the output at those energies. Table|4]lists the 
dominant coolers in the various zones. For the envelope, we use 
the temperatures listed in Table [3] 

We now discuss the individual zones in more detail to un- 
derstand the emission lines that are produced. In Sect. 14.51 we 
discuss how the energy deposition in the various zones changes 
if there is a non-local positron deposition. 



4.2.1. Fe/He zone 

The high electron fraction in the Fe/He zone (~0. 17) implies that 
a large part (~59%) of the deposited energy goes into heating, re- 
emerging mainly in the [Fe II] 26 /vm line (Table 5). The reason 
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Table 5. Energy deposition, temperature, and electron fraction in the various zones in the standard model with M( 44 Ti) = 1.5 x 
1(T 4 M . 



Zone 


rc/nc 


Jl/J 




W/ IN c/ IVlg 


\ t \ 


I 1^ 


u 

1 1 


Deposition [ergs -1 ] 


1.6 x lG ib 


9.1 x 10 34 


3.6 x 10 34 


5.4 x 10 34 


1.7 x 10 34 


2.6 x 10 34 " 


1.5 x 10 35fl 


Fractional deposition 


80% 


4.7% 


1.9% 


2.8% 


0.88% 


1.4%" 


8.0%" 


Deposition external [ergs -1 ] 


1.7 x 10 34 


7.4 x 10 33 


1.1 x 10 33 


1.1 x 10 34 


2.4 x 10 33 


7.7 x 10 33fl 


1.3 x 10 34fl 


Heating 


0.59 


0.31 


0.23 


0.19 


0.22 


0.20" 


0.14" 


Excitation 


0.16 


0.27 


0.21 


0.21 


0.18 


0.23* 


0.41* 


Ionization 


0.25 


0.42 


0.56 


0.60 


0.60 


0.57* 


0.45* 


Temperature [K] 


170 


140 


140 


70 


70 


140" 


100* 


Electron fraction 


0.17 


0.020 


2.9 x 10~ 3 


9.0 x 10~ 4 


5.6 x 10~ 4 


2.1 x 10~ 4 * 


1.5 x 10"** 



Notes. 'Deposition external' is the amount of absorbed radiation originating in other zones. 
{a) All zones. (b) Core zone. 



Table 4. Dominant cooling transitions for each zone. 



Zone 


Cooler 


Fraction 


Fe/He 


[Fe II] 25.99 fim 


94% 




[Fe II] 14.98 urn* 


5% 




[Fe I] 24.04 //m 


1.4% 


Si/S 


[Si I] 68.47 yum 


63% 




[Si I] 44.81 fim 


22% 




[Si I] 129.68 urn 


15% 


O/Si/S 


[O I] 63.19 yum 


35% 




[Si I] 68.47 yum 


32% 




[O I] 44.06 fim* 


17% 


O/Ne/Mg 


[O I] 44.06 yum 


33% 




[O I] 63.19 yum 


29% 




[Si I] 68.47 yum 


16% 


O/C 


[O I] 63.19 yum 


41% 




[O I] 44.06 yum* 


40% 




[Si I] 129.68 //m 


7% 


He (core) 


[Fell] 25.99 //m 


46% 




[Si I] 68.47 yum 


23% 




[Sill] 34.81 //m 


16% 


H (core) 


[Si II] 34.81 //m 


59% 




[O I] 63.19 yum 


19% 




[Fe II] 25.99 //m 


15% 



Notes. Transitions that are strictly forbidden (and therefore give no line 
emission) are marked with an asterisk. 



that this line dominates over [Fe I] 24 yum i s that it has a much 
larger collision strength, ~6 v ersus ~0.02 dPelan & Berringtonl 
fl997tlZhang & Pradhanll 19951) . Because the Fe/He zone absorbs 
~80% of the deposited energy (assuming on-the-spot positron 
absorption), and ~60% of this goes into heating, ~50% of the 
total deposited energy or ~ 9 x 10 35 erg s -1 (or a factor two or 
so less considering the dust) should emerge in the [Fe II] 26 /im 
line at this epoch. We discuss observations of this line in Sect[5] 
About 25% of the energy goes into ionizations, with ~14% 
going to He I, ~7% to Fe I and ~4% to Fe II. Iron reaches an 
ionization degree of ~40%, while helium reaches ~6%. Fe III is 
largely neutralized by charge transfer reactions, mainly with Fe 
I. Nickel, which makes up ~4% of the zone mass, is almost fully 
ionized owing to charge transfer with Fe II. He I, Fe I, and Ni 
I therefore emit strong recombination lines. The recombination 
line spectrum of Fe I is rather smooth with significant emissivity 
in the 2000-6000 A range; this constitutes a major part of the 
emitted flux in this range (Sect. 14.3b . The He I recombination 
emission is mainly in the form of the two-photon continuum plus 
He I A626. This flux is mostly continuum-absorbed by other ele- 



ments and serves to additionally boost the ionization of iron and 
other metals (see also KF98 b). In the optical, the strongest re- 
combination line is the 5876 A line. Helium also emits IR lines 
at 1.083 //m and 2.058 //m. The 1.083 yum line is mainly excited 
by recombinations, while the 2.058 yum line arises mainly as a 
result of non-thermal excitations. 

The remaining ~16% of the deposited energy excites reso- 
nance lines from the ground states of Fe I (~9%), Fe II (~3%) 
and He I (~3%). For the Fe I transitions, we use the approximate 
Bethe rates, whereas for Fe II and He I we have more accurate 
cross sections. 

Based on m odeling of the Fe II lines during the first 800 days 
dLi et all 1993b . the density of the Fe/He zone may be lower than 
the one we use here, but not higher. Lowering the density of this 
zone by a factor of two increases the heating fraction by about 
five percentage units, at the expense of excitations and ioniza- 
tions. Assuming that the UV/optical/NIR flux scales in propor- 
tion to the excitation plus ionization energy, the UV/optical/NIR 
emission from the Fe/He zone then falls by ~12% by such a 
maneuver, and the total UV/optical/NIR emission by ~5-10%. 
Using y core = 2500 km s _1 instead of 2000 km s (which lowers 
all densities by a factor of ~2) lowers the total ejecta emission in 
the 2000-8000 A range by -10%. 

4.2.2. Si/S zone 

The higher density and the smaller amount of 44 Ti (~5% of the 
total) present in this zone compared to the Fe/He zone implies a 
lower electron fraction (~0.02). A smaller part of the deposited 
energy therefore goes to heating (31%). Of the rest, 42% goes 
into ionizations, mainly of the abundant Si I and S I. Still, these 
elements stay almost fully neutral because of rapid charge trans- 
fer with elements such as Fe I, Mg I and Ti I. These recombine 
by further charge transfer reactions, and it is eventually calcium 
that is the dominantly ionized species. This means that quite a 
strong Ca I recombination spectrum is emitted from this zone. 
In Sect. 14.71 we show that this holds true even without charge 
transfer, because calcium is also strongly photoionized. 

About 27% of the energy goes to non-thermal excitations, 
mainly for transitions in Si I and S I. 

4.2.3. O/Si/S zone 

This is the smallest oxygen zone (0.16 M ) but contains some 
44 Ti (~2% of the total). The largest part of the energy (56%) 
goes into ionizations, which through charge transfer reactions 
eventually leave iron and aluminium as the dominantly ionized 
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elements. The non-thermal excitations occur in similar amounts 
for O I, Si I and Mg I. 



4.2.4. O/Ne/Mg zone 

The O/Ne/Mg-zone is by far the most massive oxygen zone 
(1.89 M ), but contain no 44 Ti. The total energy deposition here 
is 2.8% of the total. About half of the energy goes into ioniza- 
tions of oxygen. A series of charge transfer reactions neutralize 
the oxygen ions and transfer the ionization to sodium, which be- 
comes the most abundant ion. This zone therefore emits strong 
Na I recombination lines, including the doublet at 5890, 5896 
A as well as at 5683, 5688 A, 6154 A and 8183, 8194 A. 

In our model, the dominant charge transfer reactant with O 
II is C I, where a reaction to the excited 2p('D) state in O I is 
close to resonance with an energy de fect of only 0.4 eV. With 
the rec ipe for estimating rates given in lPequignot & Aldrovandil 
(QUI, we estimate the reaction rate to be 1 x 10 cm s , which 
together with the relatively high carbon abundance in this zone 
(4.3% by number) makes it the dominant reaction for neutraliz- 
ing the oxygen ions. Because 2p('D) is the upper state for the 
[O I] /M6300, 6364 lines, this has the interesting effect of max- 
imizing the doublet emission efficiency by producing a photon 
for each O I ionization, which is equivalent to ~7% of the de- 
posited energy. Only ~ 10% of radiative recombinations pass 
through the 2p('D) state at these densities (giving ~0.7% of the 
deposited energy), and since other competing charge transfer re- 
actions likely do not occur to 2p('D), the conversion efficiency 
is even lower if they would dominate. 

The collisional cross sections for the [O I] doublet lines 
are moderate and only ~1.5% of the deposited energy goes 
to direct non-thermal excitation of 2p('D). This state also re- 
ceives some indirect contribution by non-thermal excitations to 
higher levels, the most important channel being via the 3d( 3 S) 
state. Collisions to this state contribute ~0.4% to 2p('D) by 
cascading. The total energy emitted in [O I] AA63Q0, 6364 is 
then -8.9% (7%+1.5%+0.4%) of the deposited energy if the 
charge transfer reaction with C I dominates, and only ~2.6% 
(0.7%+1.5%+0.4%) otherwise. 

Magnesium has an abundance (by number) that is ~25 times 
lower than that of oxygen. But because the Mg I] /14571 tran- 
sition has a large collisional cross-section, it receives about the 
same number of non-thermal excitations as the [O I] doublet. 
About 3.7% of the deposited energy goes to ionizing Mg I, but 
Mg II is neutralized by charge transfer with Ni I, which in turn 
ionizes Na I. 

Sodium is five times less abundant still, but Na I has an even 
larger cross-section to its first excited state, which also receives 
~1.5% of the energy through non-thermal excitations. This con- 
tributes further to the Na I /1/15890, 5896 emission, which is, 
however, dominated by recombinations. 



4.2.5. O/C zone 

This is the intermediate-mass oxygen zone (0.58 M ), which re- 
ceives 0.9% of the total energy deposition. The low abundance 
of sodium in this zone, compared to the O/Ne/Mg zone, makes 
magnesium the dominantly ionized element. A strong Mg I] 
/14571 recombination line is therefore emitted from this zone, 
although more magnesium resides in the O/Ne/Mg zone. 



4.2.6. He zone 

Because of mixing, the helium zone has both a core component 
and an envelope component (Sect. [2]). These components behave 
similarly. He II is neutralized by charge transfer, mainly with Si 
I, and ionizations eventually end up in Na II, Mg II, Si II, and Fe 
II. 

Non-thermal excitations go mainly to the He I 2p('P) state 
(~15%), resulting in the 2.058 fim line, with only small contribu- 
tions from recombinations for this line. As long as the positrons 
are locally absorbed in the Fe/He zone, all helium emission lines 
have stronger contributions from the Fe/He zone than from the 
He zone. As an example, the fraction of the total energy de- 
posited in the ejecta going into non-thermal excitations of 2p('P) 
is 0. 16% from the helium zones, and 2.3% from the Fe/He zone, 
which means a 93% contribution to the 2.058 fim lines from the 
Fe/He zone. Similar relations hold for all other helium lines as 
well. In Sect. [5] we explore this point in more detail. Because of 
the strong helium emission from the Fe/He zone, the He zone is 
not expected to leave any particular imprint on the spectrum at 
this epoch. 

4.2.7. H zone 

As with helium, hydrogen is present both in the envelope and 
as a mixed-in core component (Sect. |2j. The envelope zones are 
powered mainly by freeze-out recombinations. Here, hydrogen 
emits a low-temperature case B recombination spectrum. The 
freeze-out emission from the envelope is generally stronger than 
the emission from the steady-state core, and it provides the bulk 
of the hydrogen lines and continua in the model. The core com- 
ponent provides the peaks of the hydrogen lines, but also serves 
to emit ionizing radiation into the other core zones. 

In the core H zone, steady-state still prevails and instanta- 
neous gamma-ray deposition determines the output. Ionization 
of H I accounts for ~36% of the deposited energy. Because we 
are close to case B, most recombinations pass through n = 2. 
About 28% of the deposited energy goes into non-thermal ex- 
citations of the H I 2p state. Both in the core and in the (inner) 
envelope, 2s and 2p are still populated according to their statis- 
tical weights, i.e. as 1:3. In the core component, the Lya optical 
depth is ~ 4x 10 8 , which gives a ratio of two-photon to Lya emis- 
sivity of ~5. The two-photon continuum therefore dominates the 
decay of the n — 2 state. The two-photon continuum can affect 
the spectrum in many ways. It can boost the ionization of metals 
with low ionization treshholds, both in the H zones themselves 
and in the other zones. It can also provide an important input 
into optically thick lines and affect the UV/optical/NIR spectrum 
by fluorescence. In wavelength windows where no or few lines 
exist, the emission may escape and contribute to the UV contin- 
uum. In the envelope, the Lya optical depth is lower and Lya 
dominates over the two-photon continuum. 

Ha from the core is powered in similar amounts by non- 
thermal excitations and recombinations (~4% of the deposited 
energy each). About 70% of reco mbinations lea d to Ha emis- 
sion in a Case B scenario at 100 K (lMartinlll988l) . 

4.3. The spectrum 

The model spectrum, for M^Ti) = 1.5 x 10~ 4 M and r d = 1, 
is shown together with the observed spectrum in Figs.|3]|5] These 
figures contain five panels. Panel (a) shows the model spectrum 
(black) overlaid on the observed spectrum (red). The observed 
spectrum was dereddened and corrected for redshift (Sect. |3]). It 
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Fig. 3. Model spectrum at 2875 days for M( 44 Ti) = 1 .5 x 1(T 4 M and r d = 1, together with the observed spectrum. The 1000-5,000 
A range, (a) The observed spectrum (red) versus the model spectrum (black). The observed spectrum was dereddened, corrected 
for redshift, and smoothed below 2200 A and in the 7500-8500 A range, (b) The contributions by various selected elements to the 
flux, (c) The part of the emerging flux that comes from scattering/fluorescence (red) and the part that is direct emission (blue), (d) 
The zones in which the photons had their last interaction; this is the zone of emission if they did not scatter and the zone where they 
last experienced scattering or fluorescence if they did. (e) The zones from where the photons were emitted before any scattering or 
fluorescence occurred. The color codings for panels 3-5 are on the right-hand side. Lya extends to to a peak flux of 1.4 x 10~ 14 erg 
s -1 A -1 cnr 2 . See also Sect.FOin the text. 



was also smoothed below 2200 A and in the 7500-8500 A range. 
Panel (b) identifies the contribution by individual elements to 
the spectrum. These are the elements that emitted the photon if 
it escaped directly, and the element that produced the last scatter- 
ing/fluorescence if it did not. Panel (c) shows the part of the flux 
that was processed by scattering and/or fluorescence in red, and 
the part that comes from direct emission in blue. Panel (d) shows 
in which zones the photons were last processed. This is the zone 
of emission for photons that escaped directly, and the zone of 
the last scattering/fluorescence for those which did not. Finally, 
panel (e) shows from which zones the photons were originally 
emitted (which is the same as in panel (d) for the photons that 
escaped directly). 

We see that the model reproduces most of the features in the 
observed spectrum, which bolsters our confidence in the funda- 
mental assumptions in the model and to the code. The main ex- 
ceptions occur in the 3000-3500 A range, which we discuss be- 
low. Apart from this range, all major observed lines are present 
in the model spectrum, and, equally important, the model does 
not produce any strong lines that are not observed. 

4.3.1. Scattering and fluorescence 

Although the line blocking effect in the photospheric phase has 
been extensively modeled before (e.g. [K arp et alJ|l977t iLucvl 
Il987t IWagoner et ail Il99lt iMazzali & Lucvlll993l) . an analy- 
sis of the ac t ual flu orescence in the nebular phase is rarer. 
iLi & McCravl d 19961) study the process at 300 and 800 days, us- 



ing a He I two-photon continuum that scatters in a single-density 
hydrogen envelope. They conclude that the process transfers a 
significant fraction of the UV emission into a quasi-continuum 
in the optical and NIR up to at least 800 days. The importance of 
the effect was found to be sensitive to the adopted temperature, 
which determines the population of the lower levels, and thereby 
the number of optically thick lines. 

That especially resonance lines can be optically thick even 
at eight years can be seen with a simple estimate. The Sobolev 
optical depth for a line from the ground state can be written as 

T - 3 x 10 3 /abs ( ( JL) (—2—) , (39) 
\5000A/U0- 5 /\10 6 cm- 3 r 

where X is the number fraction of the element and n is the to- 
tal number density, which is in the range of 10 4 - 10 7 cm -3 in 
the core at this epoch. Obviously, even elements with low abun- 
dances can still be optically thick in their ground-state resonance 
lines. 

By looking at the panels labeled (b) in Figs. 05] we see that 
scattering and fluorescence are indeed important for the spectral 
formation at eight years. The fraction of the emerging flux that 
has been processed by line transfer (red) is -60% in the UV 
(<4000 A), -30% in the optical (4000-7500 A) and -30% in the 
NIR (7500-21,000 A). Running the model with the line transfer 
switched off resulted in a -30% increase in the UV flux, a -10% 
decrease in the optical flux, and a -30% decrease in the NIR 
flux. 
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Obs. 
Model 




1 H-env 
He-env 
H-core 
He-core 
0/0 

O/Ne/Mg 
O/Si/S 
■ Si/S 
I Fe/He 
I H-env 
He-env 
H-core 
He-core 
0/0 

O/Ne/Mg 
IO/Si/S 
J Si/S 
I Fe/He 



Fig. 4. Same as Fig.|3]but for the 5000-9,000 A range. 






Unscatt. 




Scatt. 



Fig. 5. Same as Fig.[3]but for the 0.9-2.1 fim range. 



Few lines are unaffected by the line transfer. Those with the 
highest fraction of direct emission in the flux are Mg I] /14571, 
HjS, Fe I] ,15012, [Fe I] A5697, the [O I]/Fe I complex at -6300- 
6400 A, Ha+Ca 1 ,16573, He I 1.083 fim, [Si I] 1.60, 1.64 fim, 
He 1 2.058 fim, and a few others in the NIR part of the spectrum. 



the UV, has its origin in the hydrogen envelope. This is mainly 
freeze-out energy that is emitted as Lya and two-photon con- 
tinuum, and is then processed by scattering and fluorescence. 
Below ~2500 A essentially the whole spectrum originates in the 
envelope. 

By inspecting the bottom-most panels in Figs. [3]|5] we see An alternative way of illustrating the fluorescence process is 
that a substantial fraction of the emerging flux, especially in to plot the escape wavelengths versus the emitted wavelengths, 
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7000 




0.5 1 1.5 2 



Fig. 6. The fluorescence process illustrated by plotting emerg- 
ing versus emitted wavelengths (A out and A[ n , respectively) for 
all escaping photons that were absorbed by a line at least once. 
The intensities are normalized in each /l out -bin (AA/A = 0.03), so 
that the dominant Am for each A oul can be seen as the yellow/red 
regions. Contours are at 0.1, 1, 10, 50 and 90%. Very little flu- 
orescence occurs for A\ n > 7000 A, which is therefore excluded 
from the plot. 



which we did in Fig. [6] for all photons that have been scattered. 
The transformation of high-energy photons to low-energy ones 
can clearly be seen here, and it is clear that most of the fluo- 
rescent output originates as UV emission (A m < 4000 A). The 
low occupation numbers of excited states makes any inverse flu- 
orescence, i.e. conversion to shorter wavelengths, unimportant at 
this epoch, and the line transfer produces a systematic reddening 
of the spectrum. The contribution by many input wavelengths to 
each output wavelength shows that the fluoresce nce process is 
insens itive to the input spectrum, as also found bv lLi & McCravl 
(119961) . Some of the more distinct input lines are Lya and Fe 
I/Fe II emission lines at -2500 A, -2900 A an d -390 A. Lya 
powers the emission in Mg II 2795, 2802 (Sect. l4~3l2l ). 

4.3.2. Line identifications 

We now go through the spectrum, identifying and commenting 
on the most prominent lines. We use the same line peak wave- 
lengths as in C97. 

UV As panel (e) in Fig.[3]illustrates, most of the emerging far- 
UV flux (below ~ 2400 A) has its origin as emission from the 
envelope. Panel (b) shows that most of this flux has been re- 
processed by scattering and fluorescence. Longward of -2400 
A the core begins to make significant contributions. The feature 
at 2400 A is mainly scattering in the Fe II UV 2 multiplet, with 
similar contributions from the envelope and the iron zone. The 
feature is blended with an Fe I emission line on the red side. The 
feature at 2600 A is also a blend of Fe II scattering (UV 1) and 



Fe I emission. The strong line at 2513 A appears to be scattering 
in Si I UV 1, complemented by Fe I. 

The feature at -2790 A is scattering by Mg II ^2795, 2802 
in the envelope. We investigated the formation of this line and 
found it to be dominantly powered by Lya, which is absorbed 
in the Mg II 1240 A line, which then produces fluorescence at 
2795, 2802 A, 2798 A, 2928, 2936 A, 9218 A, and 10,914 A. 
There is some flux missing on the blue side of the Mg II feature, 
but we find no Fe II contributions that could improve this, as 
suggested in C97. 

Next to the Mg II feature follows emission/scattering in the 
Mg I /12852 line. This feature has a broad component that comes 
from scattering in the envelope (the line is optically thick out to 
-6000 kms -1 even though magnesium is -99% ionized in the 
envelope), and a narrow component that is dominated by emis- 
sion from the O/Ne/Mg zone. On the red side there is some con- 
tribution from iron and other elements. 

The 3000-3500 A range is not well reproduced by the model, 
and no definitive identifications can be made. Panels (d) and (e) 
in Fig. [3] show that most of the flux here originates in the Fe/He 
zone. The model produces strong lines at 2986 A (Fe I), 3090 
A (Fe I and Al I), 3160 A (Ti II scattering), 3270 A (Fe II), 3310 
A (Ni I scattering), 3400 A, 3440 A and 3490 A (emission and 
scattering mainly from iron, cobalt, and nickel). The overproduc- 
tion of several lines in this region could conceivably be due to 
the omission of trace elements that could provide additional line 
absorptions. We surveyed the line lists in the 3000-3500 A range 
of all omitted elements with Z < 30, however, and there are no 
lines from these elements that could be important scatterers at 
this epoch. The discrepancy could also be caused by the lack of 
accurate collision cross-sections for the non-thermal excitations 
in Fe I (we use the Bethe approximation). Although several of 
the lines are produced by scattering/fluorescence, the origin is 
often an Fe I emission line. As we see in Sect. 14.71 charge trans- 
fer can also be responsible by introducing uncertainties in the 
ionization balance. Finally, it is possible that the dust extinction 
is wavelength-dependent and should be higher in the UV. This 
would affect the UV lines coming from the core more than the 
ones coming from the envelope. 

We identify the feature at 3650 A as -50% Balmer contin- 
uum and -50% scattering by mainly Fe I and Cr I in the iron 
core. The 'contamination' of the Balmer continuum makes the 
attempts to determine of the hydrogen zone temperature in W96 
and C97 unsuitable, since they are based on an assumption of 
pure Balmer emission. The temperatures derived there are also 
significantly higher than in the ones calculated in our model 
(-100 K). Taking the contaminations of both the Balmer con- 
tinuum and Ha into account (see below), the observed ratio be- 
tween the Balmer continuum and Ha is 0.25, which is close to 
the theoretical ratio of 0. 32 obta i ned by using the Case B recom- 
bination rates at 80 K in lMartinl (1 1988r> . 

The strong observed line at 3730 A is reproduced by the 
model and appears to be scattering and emission from Fe I. We 
find the [O II] A3121 contribution to be negligible, agreeing with 
C97 that there is no need for this line, a possibility raised in W96. 

Optical The 4000-4500 A range is a blend of Fe I, Ca I, H I, and 
others. The strongest line in this range is at 4223 A, which the 
model identifies as Ca I /14226 emission from the Si/S zone. As 
we discussed in Sect. 14.2.21 Ca I is likely to emit most of the ion- 
ization energy in this zone, a conclusion which is supported by 
the presence of the 4226 A line in the observed spectrum. This 
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line is sensitive to density; A/3 has a value of ~600 s for the 
density we use here (a number density ten times lower than in 
the O/Ne/Mg zone), which dominates the rate of A tot = 140 s _1 
for other radiative branching paths. But if the Si/S zone would 
have the same density as the oxygen zones, A/3 would be only 
~10 s _1 (calcium is also less ionized at higher density and so 
the scaling is not proportional), and the other deexcitation paths 
would then dominate. The presence of this line in the spectrum, 
at this epoch, therefore requires a number density in the Si/S 
zone at least ~10 times lower than the one of the oxygen zones. 
Because the Si/S zone contains some (~5%) of the 56 Ni, accord- 
ing to the explosion model we use, this is consistent with the 
physical picture that radioactive decays cause expansions in the 
early phase of the su pernova, as was found to have happened for 
the Fe/He clumps bv lLi et alj d!993l) . 

The rest of the 4200-4500 A 'plateau' is mainly emission 
and scattering from the Fe/He zone, complemented with Hy at 
4343 A. The Fe II contribution, which is suggested to dominate 
the features at 4223 A, 4339 A, and 4453 A in C97, is only at 
the ~10% level in our model. 

Mg I] /14571 is one of the lines with the smallest contribution 
from scattering/fluorescence. The model underproduces the flux 
by a factor of ~2-3, possibly as a result of the uncertain ioniza- 
tion balance in the O/Ne/Mg zone (Sect. 14.2.41 ). Although most 
of the synthesized magnesium resides in this zone, its contribu- 
tion to the Mg I] /14571 line is moderate because Mg II is largely 
neutralized by charge transfer with nickel and sodium, leaving 
only non-thermal excitations to power the line. It is possible that 
we use rates that are too high for some of these charge trans- 
fer reactions, and that Mg II instead recombines radiatively, in 
which case the Mg I] /14571 line would become stronger. Our 
model gives emission of similar magnitudes from the O/Ne/Mg 
zone, the O/C zone, and the He and H zones. 

The feature at 4850 A is -2/3 H/3 and ~l/3 Fe I in the model. 
Removing the Fe I contamination gives a Balmer decrement of 
~5, which agrees well with a Case B recombination scenario at 
-100 K. 

The 5000-5700 A range is a quasi-continuum dominated by 
Fe I lines from the iron core, with about 3/4 of the flux being 
unscattered. The range is well reproduced by the model, apart 
from an overproduction of the Fe I line at 5007 A (the narrow 
observed line does not come from the ejecta, but from the north- 
ern loop). 

The line at 5900 A is emission and scattering in Na I 5896 
A. The scattered part is predominantly from emission in Na I 
5890 A, but also from He I 5876 A emitted in the Fe/He zone. 
Panel (e) in Fig. [4] shows that most of the energy originates in 
the O/Ne/Mg zone, which means that the helium contribution 
is minor. Na I emits strongly from the O/Ne/Mg zone because 
of a series of charge transfer reactions and a high photoioniza- 
tion rate that make it the dominant ion in this zone (Sect. l4~2~4b . 
Apart from the 5890, 5896 A doublet, the Na I recombination 
spectrum also includes lines at 5690 A, 6154 A, 8190 A and 
1.14 fim. There are indeed observed lines at ~5690 and ~6154 
A, but blending with Fe I and Ca I makes identifications uncer- 
tain. At 8190 A the spectrum is too noisy to clearly detect a line, 
and in addition there is also contamination from other lines in 
the model. There is a clear dete ction at 1.14 /im in the day 2112 
spectrum in lFassia e t al. (2002) (not shown here), but all in all it 
is difficult to tell whether a significant Na I recombination spec- 
trum is emitted by the ejecta or not. The model overproduces the 
Na I /1/15890, 5896 doublet, and it is possible that the ionizations 
instead go to magnesium in the O/Ne/Mg zone, especially be- 



cause the Mg I] /14571 line is underproduced. However, in Sect. 
14.71 we show that the Na I doublet is overproduced also if the 
charge transfer reactions are switched off. 

[O I] AA630Q, 6364 is contaminated, which the deviation of 
the observed ratio from the expected 3:1 value verifies. In the 
model, most of the contaminating flux is from Fe I. The complex 
is well reproduced by the model, although the good fit relies on 
the assumed high rate for the O II + C I O I(2p 'D) + C II 
charge transfer reaction, as discussed in Sect. 14.2.41 If this re- 
action is not the dominant channel for O II recombinations, the 
[O I] /L16300, 6364 emission would decrease by a factor of ~4. 
The uncertainties regarding this reaction combined with the Fe 
I contamination make the oxygen doublet unsuitable for draw- 
ing any strong conclusions at this epoch, including attempts to 
determine the oxygen mass. 

Interestingly, Ha is significantly contaminated by Ca I 6572 
A in the model. This line originates, like the Ca I 4226 A line, 
mainly in the Si/S zone, but also in the O/Si/S zone. Being the 
first transition from the ground state in Ca I, the line has a high 
effective recombination rate and is about twice as strong as the 
Ha emission from the core H zone. Note that this component is 
needed to satisfactorily reproduce the 'Hff' profile, which would 
otherwise be too flat-topped. For the high-velocity part of Ha, 
the good fit of the model is an important verification that our 
treatment of the freeze-out in the envelope is accurate. 

The [Ca II] AA729 1,7323 doublet is well reproduced and re- 
sults from pumping in the H and K lines. Fig. [4]i-e show that 
the fluorescence is produced by both synthesized and primordial 
calcium. The H, K lines are optically thick in all core zones, and 
in the envelope out to ~7000 kms~'. Fig. [6] shows that most of 
the photons start as emission close to the H and K lines. The flux 
in the Ca II /L18600 IR triplet is similar to the flux in the doublet, 
as it must be for a H, K pumping scheme. 

Infrared In the infrared, some of the observed lines are domi- 
nated by narrow circumstellar components seen at low resolu- 
tion. We i dentify the 1.083 /mi line, observed in the day 2112 
spectrum (Fass ia et al.l2002l) . with He I emission from the Fe/He 
zone. The broad c omponent has an observed flux (Table 2 in 
iFassia et all d2002l) ) similar to the model flux. Pay Al. 0938 /urn 
also has a strong narrow component. The model shows the ejecta 
emission here to be dominated by [Si I] 1 .099 jum. The model 
flux is 1.5 times the the observed broad component at 21 12 days. 

Further, we identify 1.14 fim with Na I, 1.20 jum with Si I, 
1.25 //m with [Fe II], and 1.28 fim with Pa/3. The strong 1.44 /urn 
line is [Fe I] emission from the core, with a model flux 1.8 times 
the observed value. Also 1 .50 fim and 1.53 fim are Fe I emission. 
1.60 and 1.64 //m are dominantly [Si I] emission from the Si/S 
zone, and 2.058 jum is He I from the iron core. Th e 2.058 ium 
line is overproduced in the model. We also refer to iKjaer et all 
(1201 Oh . where we discuss the IR spectrum at an age of 19 years, 
with most line identification being the same as at this epoch. 

Summary of the spectrum In general, we find good agreement 
with the line identifications in W96 and C97. The most impor- 
tant difference is that we find many lines to be from Fe I rather 
than Fe II. This is what to expect since Fe I emits not only a re- 
combination line spectrum, but also receives more non-thermal 
excitations because Fe I atoms outnumber Fe II ions by a factor 
of -1.5 (see Sect. 14.2. it . In our model, -30% of the flux be- 
tween 2000 and 8000 A is from Fe I, and only ~9% is from Fe II 
(including both emission and scattering). We do not, however, 
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include ionizations to excited states in Fe II, which causes some 
underestimate of the Fe II emission. 

4.4. The 44 77 mass 

The character of the spectrum, i.e. the line ratios and the ratio 
of lines to continuum, is not very sensitive to the absolute level 
of the radioactive powering. The character is partly determined 
by the relative amounts of energy that go into excitations and 
ionizations, which is insensitive to the total deposition. It is also 
influenced by the ionization balance of the dominant elements in 
each zone, which is likewise insensitive to the total deposition; 
the ionization fraction of the dominant element approximately 
scales as the square root of the deposition. Furthermore, there is 
a negative feedback present because a higher electron fraction 
reduces the fraction that goes into ionizations (and excitations) 
at the expense of the heating fraction. 

The 44 Ti mass therefore mainly determines the absolute flux 
level of the spectrum, rather than its form. Unfortunately, it is 
difficult to distinguish the 44 Ti mass from the uncertain dust 
properties. We therefore calculated the integrated flux in the 
2000-8000 A range for different combinations of M( 44 Ti) and 
the dust optical depth Td (Sect. l2~6l l. We did not include the NIR 
because of the uncertainties in the calibration of these observa- 
tions (Sect. [3). Fig.|7]shows the 44 Ti mass that reproduces the ob- 
served flux in the 2000-8000 A range as function of Td, together 
with the masses that produce the correct flux plus/minus 10% 
and 20%. Note that the emerging flux is less than proportional to 
the 44 Ti mass. The main reason is that a significant fraction of the 
flux comes from freeze-out emission rather than from instanta- 
neous 44 Ti deposition. Another reason is that a higher deposition 
increases the electron fraction, which increases the fraction of 
the flux emerging in the FIR. 

A 44 Ti mass of 1.5 x 10~ 4 M is our best match for Td = 1, 
which is our best estimate for the dust extinction from the line 
profiles. This result is only valid for the assumption of local 
positron deposition. The dust optical depth has a +30% uncer- 
tainty (Sect. [276b , which according to Fig. Q gives a ±20% un- 
certainty for the best fitting 44 Ti mass. In Sect. 14.2. TI we noted 
that the uncertainty in the density of the Fe/He zone adds an un- 
certainty of +10% for the model flux, and in Sect. [3] we saw that 
the observed flux level has a +15% uncertainty from errors in 
distance, reddening, and flux calibration. Together, these give a 
(15% 2 + 10% 2 ) 1/2 = 18% uncertainty for the ratio between mod- 
eled and observed flux. According to the contours in Fig. [7] this 
corresponds to a +30% uncertainty in the 44 Ti mass. Together 
with the dust error, the total uncertainty for the 44 Ti-mass be- 
comes ~ (20% 2 + 30% 2 ) 1/2 = 36%. The 44 Ti mass is then con- 
strained to (1.0 - 2.0) • 10~ 4 M . 

4.5. Positron trapping 

An described in Sect. 12.21 an interesting and important question 
is whether the positrons are indeed locally absorbed, as assumed 
so far, or if they leak into the other zones in the core before they 
are absorbed. We recomputed the spectrum for a model where 
we assume that the positrons are not absorbed on the spot, but in 
proportion to the electron content (free and bound) of each zone. 
This treatment is based on the Bethe stopping formula (Eq. ([TJ), 
where we neglect the small differences due to different effective 
ionization potentials. The number of electrons per zone is pro- 
portional to MZ/A, where M is the zone mass and Z and A are 
the average nuclear charge and atomic weight, respectively. 




.5 



Fig. 7. Relationship between the Ti-mass and the dust optical 
depth tj required to match the total luminosity in the 2000-8000 
A range (solid black line). Also plotted are the 10% and 20% 
contours of | (L mo dei - ^obs) l/^obs (dashed lines). 



For this alternative model, the Ti mass needed to give the 
correct flux in the 2000-8000 A range for r d = 1 is only 8 x 
10~ 5 M , which we therefore used. The reason for the lower 44 Ti 
mass is that a smaller fraction of the positron energy is now emit- 
ted by far-IR cooling lines (mainly [Fe II] 26 /mi) by the low- 
density Fe/He zone, and that the oxygen and hydrogen zones 
have some strong lines in the 2000-8000 A range. Consequently, 
less 44 Ti is needed to get the same UV/optical/NIR luminosity. 

In this model, the zones obtain total energy depositions 
(positrons + gamma-rays) given in Table [6] where we also give 
the relative change in the deposition fraction with respect to the 
on-the-spot model. From this, we expect positron leakage to lead 
to a severe weakening of emission lines from the Fe/He zone, 
and a boost by up to a factor of ~10 for emission lines from the 
O/Ne/Mg, O/C, core He, and core H zones. The 44 Ti content of 
the silicon zones (Si/S and O/Si/S) however, is such that the de- 
position here does not change much between the two scenarios. 

Table 6. Fractional energy deposition per core zone in the 
positron-leakage model (which roughly corresponds to the frac- 
tion of the electrons that reside in each zone), and the change 
in fractional deposition relative to the on-the-spot model. About 
5% of the total deposition occurs outside the core. 



Zone 


Fractional deposition 


Relative change 


Fe/He 


1.3% 


0.017 


Si/S 


2.1% 


0.43 


O/Si/S 


2.4% 


1.2 


O/Ne/Mg 


28% 


10 


O/C 


8.7% 


10 


He(core) 


4.9% 


10 


H(core) 


48% 


11 



Fig. [8] shows the resulting model spectrum compared to the 
standard on-the-spot model (which is the same as in Figs. |3]Q. 
On the whole, the spectra are quite similar, but contain some 
important differences as well. 

The hydrogen emission from the core is boosted by a factor 
of several, as expected. Ha and the Balmer continuum are now 
overproduced. However, because the fraction of the hydrogen 
that we mix into the core (25%) is not tightly constrained, there 
is some room to improve the H lines by reducing this in-mixing. 
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That would, however, increase the deposition in the other core 
zones and strengthen the lines from them. 

The [O I] /t/16300, 6364 doublet becomes much too strong, 
but this result is sensitive to the uncertain OII + CI— >OI 
(2p'D) + C II charge transfer rate used, as discussed in Sect. 
14.2.41 If this reaction is switched off, the agreement is actually 
better in this model than in the on-the-spot model. The Mg I] 
/14571 line is mainly emitted from the oxygen zones, and be- 
comes about as much overproduced as it was underproduced be- 
fore. The Mg I /12852 line, however, now becomes much too 
strong (it is truncated in the figure). This line is powered by 
non-thermal excitations, and is therefore independent of the un- 
certain ionization balance in the O/Ne/Mg zone (magnesium is 
dominantly neutral in all scenarios). This line should therefore 
be a more reliable indicator than the recombination-driven 4571 
A line, and all this suggests that too much energy is now being 
deposited in the oxygen zones. The Na I recombination lines at 
5890, 5896 A, and 6160 A from the O/Ne/Mg zone become over- 
produced as well. The 5890, 5896 A doublet, which was already 
too strong without leakage, is now a factor ~4 too strong. 

For the rest of the spectrum, the two models produce sur- 
prisingly similar spectra. Because fluorescence makes impor- 
tant contributions at almost all wavelengths, turning off the iron 
emission lines from the Fe/He zone typically lowers the line 
fluxes by no more than a factor of about two. An example of 
this situation is the 5000-5500 A range. The iron lines with the 
smallest contamination by fluorescence in the UV/optical are the 
lines at 3270 A, 3490 A, 5012 A, and 5159 A, and 5697 A. All 
of these lines seem somewhat better reproduced in the leakage 
model. 

Because the flux calibration in the NIR is uncertain (Sect. 
[3]), we refrain from comparing the models in this range. All in 
all, we think the on-the-spot model produces a spectrum that 
better agrees with the observations, although the uncertainties 
mentioned above prohibit a definitive ruling out of the leakage 
scenario. The average x 1 -value (per spectral bin) is about twice 
as high in the leakage model as in the on-the-spot model. 

4.6. Fragmentation level 

The fragmentation level, N c \, is a proxy for the hydrodynami- 
cal mixing. This parameter determines the length of the paths 
that the photons travel in their parent zone before they have the 
chance to enter the other core zones. A high fragmentation level, 
implying quick exposure to other zones, will lead to a higher 
degree of radiative energy exchange between the zones. If one 
were able to constrain this parameter it would be very useful for 
evaluating multi-dimensional explosion models. 

To check the sensitivity of the spectrum to this parameter, 
we compared spectra for models using N c \ = 10 and N c i = 10 5 . 
The results were very similar, and we therefore conclude that 
the fragmentation parameter has no significant influence on the 
spectrum at this epoch. One possible explanation for this is that 
few line optical depths can be expected to be close enough to 
unity that an abundance change (by switching zones) between a 
high value in the zones where the element is synthesized to the 
(low) primordial values where it is not, will make any difference. 
Another could be that few of the line absorptions occur before 
the photon has exited its parent zone, even in the low fragmen- 
tation limit. Finally, a significant part of the line transfer occurs 
in the envelope, where there is no mixing in the model. 
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Fig. 8. Standard model (assuming on-the-spot positron absorp- 
tion) (upper, black) compared to a model assuming positron de- 
position in proportion to the electron content of each core zone 
(lower, black). The latter uses a 40% smaller 44 Ti mass (in order 
to produce the correct total flux in the 2000-8000 A range). The 
Mg I 2852 A line continues up to ~ 3.1 x 10~ 15 in the bottom 
panel. The observed spectrum in red, same as in Fig. [3] 



4. 7. Effects of charge transfer 

Finally, the uncertainty in many of the charge transfer rates 
prompted us to examine the influence of these reactions on the 
spectrum. It is also interesting to see how important these reac- 
tions are for supernova spectra in general. We recomputed our 
standard model (M(^Ti) = 1.5 x 10~ 4 M , r d = 1, and on- 
the-spot positron deposition) with all charge transfer reactions 
switched off. The resulting spectrum is compared to the standard 
one in Fig. [9] 

The total flux in the 2000-8000 A range changes by only 
a few percent, so charge transfer should not affect our conclu- 
sions regarding the 44 Ti mass. Several lines differ significantly 
between the two models though. The [O I] /L16300, 6364 dou- 
blet is weakened for reasons discussed in Sect. 14.2.41 and the 
Fe I lines now dominate the doublet. At the same time, the O I 
/17775 recombination line, which is not observed, emerges in the 
model. This suggests that the oxygen ions are indeed neutralized 
by charge transfer reactions, as in our standard model. 

The important Ca I lines (4226 A and 6572 A) do not de- 
crease by much, which shows that charge transfer is not criti- 
cal for them. Ca I is significantly ionized by the internal radia- 
tion field, and so emits strong recombination lines even without 
charge transfer. 

The Na I /15896 emission is reduced somewhat, but is still 
too strong. Just as with Ca I, Na I is significantly ionized by the 
internal radiation field, and consequently strong recombination 
emission still occurs. An important conclusion is therefore that 
the emission lines from elements with low ionization thresholds 
are is not very sensitive to the charge transfer, because low ion- 
ization thresholds also correlate with high radiative ionization 
rates. 

For many lines, it is not possible to directly say why they 
change. The 'picket fence' of absorption lines changes as the 
ionization balance changes, which opens some wavelength win- 
dows and closes others. Consider, for example, the Fe I emission 
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Fig. 9. The standard model (upper, black) compared to a model 
with all charge transfer reactions switched off (lower, black). 
Observed spectrum in red. 



line at ~2990 A. This emission line is still there in the model 
without charge transfer, but now a Ni I absorption line has be- 
come optically thick that absorbs it. The UV spectrum is actu- 
ally somewhat better reproduced in the model without charge 
transfer, which may indicate that we overestimate some of our 
adopted rates. Unknown charge transfer rates clearly constitute 
one of the main obstacles to accurate supernova spectral mod- 
eling today, and more calculations and measurements of these 
would be highly desirable. 



5. Discussion 

Although most of the optical and infrared emission escapes 
freely in the nebular phase of supernovae, the UV emission 
does not. The freeze-out energy in the envelope, as well as non- 
thermal excitations and ionizations in the core, produce strong 
UV emission that is significantly absorbed by atomic lines and 
emerges by fluorescence at longer wavelengths. By this process, 
the emerging flux in the optical and NIR is enhanced by ~10% 
and ~30% compared to a purely nebular spectrum at the epoch 
we have looked at here. While the non-local line transfer can 
have a large impact on individual lines at late times, the fact that 
it has a moderate impact on the supernova photometry is impor- 
tant forjudging the accuracy of models without this effect, which 
has previously been difficult. 

While the fraction of UV photons that are absorbed de- 
creases with time, the fraction of the optical/NIR spectrum that 
come s from fluorescence may actually increase dLi & McCravl 
119961) . The reason for this is that with time the cooling lines 
move from the optical/NIR into the FIR, and non-thermal pro- 
cesses do not produce very strong optical/NIR emission. Our 
finding that fluorescence is important also at eight years illus- 
trates this point well. 

We showed that the only FUV flux that escapes comes 
from far out in the hydrogen envelope. A consequence of this 
is that the ejecta appear larger in the U V than in the optical 
dJakobsen et al.lll994l: IWang et al1ll996l) . At longer UV wave- 
lengths both the envelope and the core contribute to the spec- 
trum. 



From a mathematical perspective, it is not a very well- 
conditioned problem to attempt exact modeling of the UV lines, 
because the output depends sensitively on the exact shape of 
the 'picket fence' of absorption lines, which in turn depends on 
many uncertain parameters such as abundances, densities, and 
charge transfer rates. Considering this complexity, it is probably 
not realistic or meaningful to attempt any modeling of the UV 
lines in much more detail than we have done here. 

As discussed in Sect. 12.21 the positrons emitted in the 44 Ti 
decay will enter the neighboring oxygen, helium, and hydrogen 
zones at late times unless a magnetic field exists to trap them. 
This would have significant consequences for the FIR cooling 
lines. The model flux in the [Fe II] 26 /urn line, discussed in 
Sect. 14.2.11 corresponds to a peak flux of ~1.8 Jy at day 2875. 
ISO observations at day 3425 failed to detect this line, giving a 
upper limit for the peak flux of 0. 64 Jv ([Lundavist et al.ll200ll) . 
Spitzer observations at day 6190 dBouchet et al.l 120061) show a 
distinct broad feature (FWHM ~ 3800 kmr 1 ) at 26 //m, but 
with a peak flux of only ~0.02 Jy. The slow decay of 44 Ti should 
make the situation at 6190 days similar to the one at 2875 days, 
and it is therefore surprising that the observed flux in this line is 
about two orders of magnitude weaker than expected. Possible 
explanations are that the positrons are indeed leaking out into the 
other core zones, or that dust is cooling the Fe/He zone. In the 
positron leakage model, the flux in the 26 //m line is ~0. 16 Jy at 
2875 days (for a 44 Ti mass of 8 x 10~ 5 M and no dust correc- 
tion). As we found in Sect. 14.51 a positron leakage scenario pro- 
duces a UV/optical spectrum with several lines that are clearly 
too strong. Our favored explanation for the weak [Fe II] 26 //m 
line is therefore that dust cooling occurs in the Fe/He zone. The 
possibilities of h aving dust cooling in t he ejecta was discussed 
in KF98 a and in lLundqvist et all (1200 lh . 

In the case of magnetic confinement of the positrons, all he- 
lium lines are dominated by helium in the Fe/He zone, which 
is produced by freeze-out in the explosion. Although we do not 
identify any distinct observed lines with helium in the eight-year 
spectrum, the 19-year spectrum in dKjaer et al 1 120101) contains a 
strong and distinct 2.058 fim line, which, as we speculated there, 
probably shows that ff-rich freeze-out has indeed happened in 
the Fe/He zone; without it the local positron deposition could 
not produce such a strong line. 

Together with the 56 Ni mass of 0.069 Mn dBouchet et al. 
1991b and the 57 Ni mass of ~ 3.3 x 10 3 M m dFransson & Kozma 



1993blVarani et al.lll990UKurfess et aDl 19921) . the determination 
of the 44 Ti mass should put considerable constraints on super- 
nova explosion models. Our result of 1.5^5 - 5 x 10~ 4 M Q agrees 
with the result of (1 - 2) x 10~ 4 M found in C97, although that 
result was based on an analysis of [Fe II] lines, whose identifi- 
cations are different in our model on several instances, and on 
the [O I] /L16300, 6364 doublet, whose emission is highly sen- 
sitive to uncertain charge transfer reactions, and which is also 
contaminated by Fe I lines. It is therefore a bit of a coincidence 
that w e get similar results for the 44 Ti mass. lFransson & Kozmal 
d2002l) determine the 44 Ti mass to (0.5-2) x 10~ 4 M Q by compar- 
ing the observed B and V band HST photometry between 1800 
and 3500 days to the output from a time-dependent model with- 
out non-local radiative transfer (as in KF98 a, b). We showed 
here that this transfer, although it is important for individual 
lines, only boosts the output in the optical/NIR range by 10- 
30% at 2875 days (Sect. 14.3.11 ), which implies that the results 
in FK02 should have a corresponding accuracy. An inspection 
of Fig. 4 in the FK02 paper shows that the best-fitting 44 Ti mass 
is in the (1.5 - 2) x 10~ 4 M range, which agrees well with our 
results here; ignoring non-local radiative transfer should give a 
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somewhat over-estimated 44 Ti mass, but not by much. Our result 
also agrees with the range o f (0.8 - 2.3) x 1CT 4 M found by 
iMotizuki & Kumagail j2004t) . Their result is, however, is based 
on a comparison betwen the deposited energy and the bolometric 
luminosity at 3600 days, the latter being highly uncertain since 
most of the flux comes out in the unobservable FIR. 

The 44 Ti isotop e is a result of the a-rich freeze-out (e.g. 
IWooslev etalJ 19731) . which is a consequence of the slow rate for 
the triple-alpha reaction, which cannot keep the 4 He abundance 
at the NSE value as the temperature and density decreases. The 
high 4 He abundance in turn results in a higher abundance of 44 Ti 
than is the case in NSE. The 4 4 Ti/ 56 Ni ratio therefore sens itively 
depends on the entropy (e.g. IWooslev & Hoffman! 1 19911) . High 
temperature and/or low density increases this ratio. The explo- 
sion model we use (from WH07) has a mass ratio 44 Ti/ 56 Ni = 
3.1 x 10~ 4 , which gives a 44 Ti mass of only 2.1 x 10~ 5 M (for 
M( 56 Ni) = 0.069 M ), which is clearly ruled out according to 
our modeling. For a spherically symmetric explosion of a 20 
M Q star, WW95 fi nd an ejected ma s s of ( 1 - 2) x 10~ 5 M Q . 
On the other hand, iThielemann et al.l (1 19961) find a much larger 
mass of 1.7 x 10~ 4 M Q , closer to the value we derive here. 
The differences are partly caused by the different nuclear reac- 



-driven), while 


Thielemann et alJ 


(iTimmes et al. 


19961). The solar 



44 Ti-production 
3, averaged over all progeni- 



abundance of the decay product 44 Ca requires a 
higher than in WW95 by a factor 
tor masses, suggesting th at the yields in those models are indeed 
lo w (ITimmes et al II 19961) . 

iNagataki et al] (1 19981) . iNagatakil d2000l) and 
iMaeda & Nomotd d2003l) have investigated the effects of 
asymmetries in the explosion. These (simplified) models have 
considerably higher temperatures and entropies in the polar 
direction. Consequently, the 44 Ti production occurs mainly in 
this direction, and the total 44 Ti mass can be considerably higher 
than in a spherically symmetric model. 

Further, the rates of several of th e nuclear reactions i nvolved 
in the 44 Ti production are uncertain. iNassar etaL 1 (120061) find an 
increased rate for the important 40 Ca(o', y) 44 Ti reaction, resulting 
in a factor of two higher 44 Ti abundance. The theoretic al range of 
44 Ti m ass is therefore highly uncertain, as discussed bv lThe et alJ 
d2006l) for Cas A. 

In addition to 56 Ni, 57 Ni and 44 Ti, other radioactive isotopes 
may be present in the ejecta. In particular, 6<) Co and 22 Na may 
have sizeable abundances. These isotopes are produced by neu- 
tron and proton capture, respectively, during carbon burning, and 
have unc ertain yields that ar e sensitive to the convective pre- 
scription (ITimmes et al.lll996T) . 

The 6l) Co-isotope has a decay time of 7.605 years, and may 
have a similar mass to 44 Ti (computed as 2.5 x 10~ 5 M Q in 
WW95 and 6.5 x 10" 5 M Q in WH07). The gamma-ray luminos- 
ity for these two masses are 56% and 150% of the 44 Ti gamma- 
ray luminosity, respectively, and the electron deposition is 12% 
and 30% of the 44 Ti positron deposition, respectively. 60 Co may 
therefore make a non-neglegible contribution to the powering at 
eight years. The jfi-decays occur as electrons with energies up 
to 317 keV, with an average energy of 96.5 keV. This is signifi- 
cantly lower than the energy of 44 Ti-positrons, and the stopping 
range is therefore much shorter. A major difference is that most 
of the 60 Co electron energy should be deposited in the oxygen- 
rich zones. These account for approximately half of the total col- 
umn density of the core, and from Fig.|2]it can be seen that they 
can easily trap most of these electrons. The increased deposition 
in the oxygen-rich regions means that an appreciable 60 Co mass 



may mimic the effects of the positron leakage in the 44 Ti dom- 
inated models (Sect. 14.51 ). It will in particular increase the O I, 
Mg I, and Na I line fluxes noticeably, further complicating the 
question of where the positrons are deposited. 

22 Na has a decay time of 3.75 years, emitting a gamma-ray 
at 1.274 MeV and in 90% of the cases a positron, while 10 % 
of the cases occur by electron capture. The positrons have an 
average energy of 215.5 keV and a maximum energy of 545.4 
ke V. The 22 Na mass in the explosion model is 4.5 x 10~ 6 M Q , 
in IWooslev etall (fl989l) it is 2.0 x lQ 6 , in WW95 it is 3.0 x 
10 7 Mm and in lThielemann et all (1 19961) it is 1.3xl0~ 7 M Q . This 
gives gamma-ray luminosities of 0.5%- 17% of the 44 Ti gamma- 
ray luminosity, and positron depositions of 0.2%-5.8% of the 
44 Ti positron deposition. Based on these models it is therefore 
unlikely that 22 Na is important for the energy budget at the epoch 
studied here. 

It is, unfortunately, difficult to constrain the contribution 
from 60 Co and 22 Na to the eight-year-spectrum because of the 
uncertainties in the modeling of the oxygen zones, and the pos- 
sibility of non-local positron/electron absorption occurring. We 
assume both of these isotopes to make neglegible contributions 
to the powering at eight years. The shorter decay time scales 
of 60 Co and 22 Na, compared to 44 Ti, means that their influence 
varies with time. The modeling of the spectrum and light curve 
at other epochs should therefore have the potential to distinguish 
between the 44 Ti and the 60 Co/ 22 Na input. Unfortunately, only a 
few spectra exist at epochs later than the one studied here, and 
they are increasingly contaminated by the flux from the ring col- 
lisions. Modeling of the NIR spectrum at 19 years has shown 
that the combination of M( 44 Ti) = 1 x 10~ 4 M Q and r d = 1 
matc hes the observed NIR spectrum at this epoch reasonably 
well dKjaer et al.ll2010l) . However, this spectrum may be affected 
both by scattered light a nd X-ray input from the circumstellar 
ring (lLarsson et al.ll2.01 H) . 

Small or moderate masses of 60 Co and 22 Na are supported 
by the modeling of the SN 1987 A photometry (FK02), which 
shows that the light curves do not deviate significantly from a 
pure 44 Ti input between 1800 and 3500 days. 

6. Conclusions 

Our main conclusions are: 

- A spectral model based on the determination of the positron 
and gamma-ray degradation, the ionization, excitation and 
temperature structure of the ejecta, and the multi-line radia- 
tive transfer, reproduces most of the observed lines in the 
eight-year UV/optical/NIR spectrum of SN 1987A to good 
accuracy. We showed that many regions in the spectrum are 
produced by Fe I lines. The contribution by Fe II lines is 
small. 

- Even many years after explosion, scattering and fluorescence 
are important processes for the spectral formation in Type- II 
supernovae. At an age of eight years, ~30% of the emerging 
optical/NIR flux and -60% of the UV flux in SN 1987A is 
produced by fluorescence. Non-local line transfer decreases 
the UV flux by ~30%, increases the optical flux by ~10%, 
and increases the NIR flux by ~30%, compared to a nebular 
treatment. 

- Most of the flux in the hydrogen lines and continua, as 
well as significant parts of the UV spectrum, are at eight 
years produced by freeze-out emission from the envelope 
rather than by instantaneous 44 Ti deposition. In the far-UV 
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(< 2500 A), essentially the whole spectrum is produced in 
the envelope. 

- A scenario of positron leakage from the Fe/He clumps pro- 
duces too strong emission in several lines from the hydrogen 
and oxygen z ones. We therefore support the conclusion in 
IChugai et al.l d 1997b that the positrons remain trapped by a 
magnetic field. 

- The [Fe II] 26 jum line in Spitzer observations at day 6190 is 
much weaker than expected. Possible explanations are dust 
cooling or that positron leakage is occuring at that epoch. 

- Modeling the dust as a gray absorber with radial optical 
depth Td = 1, which is our best estimate from the line peak 
blue-shifts, the amount of 44 Ti synthesized in the explosion 
that best reproduces the total flux between 2000 and 8000 
A is 1.5 x 10 -4 M . We estimate the uncertainty in this result 
to +0.5 x 10~ 4 M . 

- One of our line identifications is the Ca I AA226 line. If this 
identification is correct, the density in the Si/S zone must be 
at least ~10 times lower than in the oxygen zones for the line 
to be as strong as observed at this epoch. Because the Si/S 
zone contains some of the 56 Ni, this adds further evidence to 
the picture that the early radioactive decays cause expansion 
in the ejecta. 
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Appendix A: Atomic data 

Table I A. 1 1 lists a summary of the model atoms used and the 
sources for the atomic data. The number of levels listed are the 
number of levels (including fine-structure) whose populations 
are solved for, followed by the total number of levels included 
for making the line list. For the population solutions, a size of 
50 levels for the neutrals and 10 levels for the ions were typi- 
cally used, with the larger number for the neutrals because these 
emit recombination spectra, whereas no ion does in any signifi- 
cance. No lines with optical depths higher than 10~ 3 have their 
lower levels above these limits at this epoch. Elements that have 
non-thermal excitations or specific recombination rates to higher 
levels, have more levels included in the population solution. Two 
examples are Fe I and Fe II, where the populations of all ~500 
levels are calculated. 

All atoms are solved for with fine-structure, which is es- 
pecially appropriate at late times. For allowed transitions that 
lack pub lished collision strengt hs, we use the van Regemorter 
formula (Ivan Regemorterlll962l) . For forbidden transitions we 
used the van Regemorter formula with an absorption oscillator 
strength of 1 x 10~ 2 g 1( . This treatment is based on the fact that 
forbidden transitions usually have collision strengths about an 
order of magnitude smaller than allowed transitions. This treat- 
ment also allows a reproduction of the typically very large colli- 
sion strengths for transitions of the type nl — > «/', since the van 
Regemorter formula contains a AE' 1 dependence. We checked 
the sensitivity to this treatment by comparing the spectrum with 
one using a constant T u j = Q.lgig,, for the forbidden lines. Only 
the Mg I] /14571 line was sensitive to the treatment. The fine- 
stru cture collision strength s in this atom is on the order of 10 5 
(e.g. lSundqvist et al. 2008), and the first treatment is better. 



If specific recombination rates are lacking, we allocate 30% 
of the total recombination rate to the ground multiplet and the 
remaining 70% in proportion to the statistical weights. If some 
specific recombination rates are available but not all, we allocate 
the difference between the total rate and the sum of the known 
rates to the states with unknown rates in proportion to their sta- 
tistical weights. If all specific recombination rates are known, but 
do not match the total rate, we scale them so that they do. We do 
not include dielectronic recombinations, which are unimportant 
at the low temperatures here. 

All ground-s tate photo-ionizatio n cross sections are from a 
routine based on lVerner et al.l (1 1996b . In addition, cross sections 
for some selected excited states are included (n=2 in H I, 2'S 
and 2 1 ? in He I, ! D, 'S and 5 S in C I, 'D and >S in O I, 3 P in 
Mg I, 2 P in Mg II, 'D in Si I, ! D in S I, 2 D in Ca II, a 5 F in Fe I 
and a 4 F in Fe II), taken from TOPBASFfl 

For information on the cross sections used in the Spencer- 
Fano routine, see KF92 and K98 a, b. An imp ortant update is the 
addition of detailed cross sections for Fe II (iRamsbottom et all 
l2005Ll2007h . and for Ca I (ISamson & Berringtonl200ilh 
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Table A.l. Summary of the model atoms used and sources for the atomic data. Col. 1: Element, Col. 2: number of fine-structure 
levels solved for and the number included in the radiative transfer, Col. 3: highest state included, Col. 4: source for level energies and 
transition A-values, Col. 5: source for the total radiative recombination rate, Col. 6: source for the specific radiative recombination 
rates, and the number of levels for which these rates are included in parenthesis, Col. 7: source for the (thermal) collision strengths, 
and the number of levels with collision strength included in parenthesis. 



Element 


Levels 


Highest state 


Energies and A-values 


Tot. rr. 


Spec. rr. 


Coll. strengths 


HI 


210/465 


n=20/n=30 


23 


18 


4 (all) 


1, 5 (all) 


He I 


29/29 


5p('P) 


23 


6 


3 (all) 


25 (all) 


CI 


90/90 


5d('D) 


23 


11 


11 (first 17) 


16 (first 6) 


CII 


7/77 


8g( 2 G) 


23 


10 


- 


16 (first 5) 


NI 


68/68 


5s( 2 P) 


23 


19 


- 


16 (first 5) 


Nil 


10/148 


2p 3 ( 5 S)/3p( 3 D) 


23 


19 


- 


16 (all) 


OI 


135/135 


8d( 3 D*) 


23 


10 


10 (first 49) 


16 (first 7) 


Oil 


16/16 


2p 4 ( 2 S) 


23 


10 


- 


16 (first 5) 


Nel 


3/3 


3s( 2 [3/2]) 


23 


19 


— 


- 


Nell 


3/3 


2p 6 ( 2 S) 


23 


19 




16 (first 1) 


Nal 


66/66 


9f( 2 F) 


23 


14 


24 (first 60) 


20 (first 2) 


Mgl 


200/371 


37d( 3 D)/80d( 3 D) 


7 


14, 19 


24 (first 18) 


8 (first 21) 


Mg II 


35/35 


7p( 2 P) 


23 


21 


- 


16 (first 2) 


All 


50/86 


10d( 2 D)/21s( 2 S) 


7 


14 


- 


C 


A1II 


10/154 


3p 2 ( 3 P)/3pd( 3 D) 


7 


_b 


- 


_C 


Si I 


200/494 


10s( 3 P)/21s('P) 


7 


n 


11 (first 49) 


c 


Sill 


77/77 


9p( 2 P) 


23 


10 


- 


16 (all) 


SI 


125/125 


8f( 3 F) 


7 


14, 19 


- 


C 


SII 


5/5 


3p 3 ( 2 P) 


23 


10 


- 


16 (all with gs) 


Arl 


3/3 


3p 5 ( 4 S) 


23 


19 


- 


c 


Aril 


2/2 


3p 5 ( 2 P) 


23 


19 


- 


16 


Cal 


198/198 


6d( 3 D) 


7 


14 


24 (first 11) 


c 


Call 


69/69 


16d( 2 D) 


23 


14 


- 


9 (all) 


Scl 


50/254 


4p( 4 F)/7p( 2 P) 


7 


a 


- 


C 


Sc II 


10/165 


3d 2 ( 3 P)/6f('H) 


7 


b 


- 




Til 


50/394 


4s(c 3 P)/4p 2 ( 1 P) 


7 


14 


- 


c 


Till 


10/212 


4s(a 2 F)/5d( 4 H) 


7 


_b 


- 


_c 


VI 


50/502 


sp(z 6 D)/4p 2 (r 2 H) 


7 


_a 


- 


_c 


VII 


10/323 


4s(a 5 F)/4f(2(ll/2)) 


7 






c 


CrI 


50/392 


sp(z 7 D)/21p(e 5 F) 


7 


a 




V 


Cr II 


10/500 


4s(a 4 D)/4f( 4 H) 


7 


_b 




c 


Mnl 


50/293 


4s(b 4 G)/4p(z 2 I) 


7 


14 




_c 


Mn II 


10/509 


4s(a 5 G)/9h( 7 H) 


7 


j> 




c 


Fel 


496/496 


4p( 1 P°)( 3 F) 


7 


10 


10 (first 356) 


2, 15 (first 120) 


Fell 


578/578 


( 7 S)4d( 8 D) 


12, 23 


10 


10 (first 381) 


22 (first 372) 


FeUI 


168/168 


4p( J D 


7 


_b 




C 


Col 


50/317 


d9(c 2 D)/5s(37°) 


7 


_a 




_c 


Coll 


10/108 


4s(b 3 F)/7d( 5 H) 


7, 17 


_b 




_c 


Nil 


50/136 


4p(F*)/4d( 2 [l/2]) 


23 


19 




_c 


Nill 


10/500 


4s( 2 D)/sp( 4 D) 


7 


19 




_c 
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